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ABSTRACT 

INVISCID  DISTURBANCE  DYNAMICS  IN  BAROTROPIC  SHEAR  FLOWS 


The  inviscid  nature  of  disturbance  evolution  in  shear  flows  is  investigated  as  an  initial-value  problem 
within  the  framework  of  nondivergent  vorticity  dynamics.  To  ensure  a  basic  understanding  of  physical  pro¬ 
cesses,  disturbance  evolution  is  first  considered  in  a  rectilinear  system  of  simple  shear.  Particular  emphasis  is 
placed  on  identifying  how  the  disturbance  evolution  depends  on  the  zonal  wavenumber  and  on  the  meridional 
structure  of  the  initial  conditions. 

Insight  acquired  from  the  rectilinear  problem  is  then  applied  to  a  bounded  Rankine  vortex.  Here,  the 
dependency  of  disturbance  evolution  on  the  azimuthal  wavenumber  is  of  special  interest.  Recent  development 
of  a  low-frequency  balance  theory  for  rapidly  rotating  vortices  has  provided  observational  evidence  that  the 
low  azimuthal  wavenumber  asymmetries,  especially  wavenumber  one,  are  dominant  in  the  near-vortex  region. 
The  results  of  this  work  provide  further  theoretical  evidence  of  an  inviscid  wavenumber  selection  mechanism 
that  preferentially  damps  the  higher  wavenumber  asymmetries. 

The  radial  structure  and  location  of  the  initial  conditions  are  found  to  be  critical  factors  in  determining 
how  rapidly  a  disturbance  is  compressed  or  elongated.  This  in  turn  controls  the  rate  of  disturbance  growth 
or  decay.  For  swirling  flows,  a  definition  of  an  effective  shear  that  accounts  for  both  the  radial  variations  in 
the  initial  conditions  as  well  as  the  radial  variation  in  the  angular  velocity  is  proposed.  Using  the  reciprocal 
of  this  effective  shear,  time  scales  for  a  disturbance  to  decay  to  half  its  initial  energy,  the  half-life  time,  are 
calculated  for  initial  conditions  and  symmetric  wind  profiles  that  are  found  in  hurricanes. 

Simple  shear  flow  and  the  bounded  Rankine  vortex  do  not  admit  discrete  modal  solutions  since  there  is 
no  mean  state  vorticity  gradient  to  support  them.  The  unbounded  Rankine  vortex  is  briefly  considered  in 
order  to  investigate  how  the  presence  of  discrete  neutral  modes  modifies  the  nonmodal  solutions  presented 
in  this  work. 


THESIS 


INVISCID  DISTURBANCE  DYNAMICS  IN  BAROTROPIC  SHEAR  FLOWS 


Submitted  by 
Gerald  B.  Smith,  II 
Department  of  Atmospheric  Science 


In  partial  fulfillment  of  the  requirements 
for  the  Degree  of  Master  of  Science 
Colorado  State  University 
Fort  Collins,  Colorado 


Summer  1994 


COLORADO  STATE  UNIVERSITY 


May  12.  1994 


WE  HEREBY  RECOMMEND  THAT  THE  THESIS  PREPARED  UNDER  OUR  SUPER¬ 
VISION  BY  GERALD  B.  SMITH,  H,  ENTITLED  INVISCID  DISTURBANCE  DYNAMICS  IN 
BAROTROPIC  SHEAR  FLOWS  BE  ACCEPTED  AS  FULFILLING  IN  PART  REQUIREMENTS 


FOR  THE  DEGREE  OF  MASTER  OF  SCIENCE. 


lOc^sJJdr 


ABSTRACT 


INVISCID  DISTURBANCE  DYNAMICS  IN  BAROTROPIC  SHEAR  FLOWS 

The  inviscid  nature  of  disturbance  evolution  in  shear  flows  is  investigated  as  an  initial-value 
problem  within  the  framework  of  nondivergent  vorticity  dynamics.  To  ensure  a  basic  understanding 
of  physical  processes,  disturbance  evolution  is  first  considered  in  a  rectilinear  system  of  simple 
shear.  Particular  emphasis  is  placed  on  identifying  how  the  disturbance  evolution  depends  on  the 
zonal  wavenumber  and  on  the  meridional  structure  of  the  initial  conditions. 

Insight  acquired  from  the  rectilinear  problem  is  then  applied  to  a  bounded  Rankine  vortex. 
Here,  the  dependency  of  disturbance  evolution  on  the  azimuthal  wavenumber  is  of  special  interest. 
Recent  development  of  a  low-frequency  balance  theory  for  rapidly  rotating  vortices  has  provided 
observational  evidence  that  the  low  azimuthal  wavenumber  asymmetries,  especially  wavenumber  one, 
are  dominant  in  the  near-vortex  region.  The  results  of  this  work  provide  further  theoretical  evidence 
of  an  inviscid  wavenumber  selection  mechanism  that  preferentially  damps  the  higher  wavenumber 
asymmetries. 

The  radial  structure  and  location  of  the  initial  conditions  are  found  to  be  critical  factors  in 
determining  how  rapidly  a  disturbance  is  compressed  or  elongated.  This  in  turn  controls  the  rate 
of  disturbance  growth  or  decay.  For  swirling  flows,  a  definition  of  an  effective  shear  that  accounts 
for  both  the  radial  variations  in  the  initial  conditions  as  well  as  the  radial  variation  in  the  angular 
velocity  is  proposed.  Using  the  reciprocal  of  this  effective  shear,  time  scales  for  a  disturbance  to 
decay  to  half  its  initial  energy,  the  half-life  time,  are  calculated  for  initial  conditions  and  symmetric 
wind  profiles  that  are  found  in  hurricanes. 
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Simple  shear  flow  and  the  bounded  Rankine  vortex  do  not  admit  discrete  modal  solutions  since 
there  is  no  mean  state  vortidty  gradient  to  support  them.  The  unbounded  Rankine  vortex  is  briefly 
considered  in  order  to  investigate  how  the  presence  of  discrete  neutral  modes  modifies  the  nonmodal 
solutions  presented  in  this  work. 
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CHAPTER  1 


INTRODUCTION 

Intense  vortices  observed  in  the  atmosphere  and  ocean  exhibit  a  high  degree  of  axisymmetry, 
although  they  are  subject  to  persistent  asymmetric  forcing.  In  hurricanes,  the  interaction  between 
the  vortex  and  environmental  asymmetries  has  implications  for  spiral  band  dynamics  (Willoughby, 
1977;  Guinn  and  Schubert,  1993),  intensification  (or  weakening)  processes  (Molinari,  1992),  and 
storm  motion  (Shapiro  and  Ooyama,  1990;  Shapiro,  1992;  and  Smith  and  Weber,  1993).  Thus, 
understanding  the  asymmetric  dynamics  of  intense  vortices  is  critical  for  understanding  the  physical 
mechanisms  controlling  vortex  evolution. 

Until  recently,  the  full  primitive  equations  have  been  the  favored  tool  for  investigating  the  three- 
dimensional  asymmetric  dynamics  of  intense  vortices.  Formulation  of  a  low-frequency  balance  theory, 
in  which  high-frequency  gravity  and  inertial  waves  are  filtered  while  retaining  the  pertinent  aspects  of 
advection,  is  complicated  for  rapidly  rotating  vortices  since  the  time  scale  of  the  tangential  advection 
is  comparable  to  that  of  the  gravity  and  inertial  waves  on  the  vortex.  Shapiro  and  Montgomery 
(1993,  hereafter  SM)  proposed  a  three-dimensional  asymmetric  balance  (AB)  theory  that  includes 
the  full  inertial  effects  of  the  rapidly  rotating  region  while  filtering  the  gravity  and  inertial  waves. 
In  the  vortex  environment,  AB  theory  reduces  to  the  quasi-geostrophic  balance  theory.  In  the 
absence  of  asymmetries,  AB  theory  reduces  to  Eliassen’s  axisymmetric  balance  model  (Eliassen, 
1951)  throughout  the  vortex. 

SM  were  able  to  separate  the  advective  processes  from  the  inertial  and  gravity  waves  by  defining 
a  local  Rossby  number  that  accounted  for  the  varying  rotation  in  an  intense  vortex.  This  local 
Rossby  number  was  given  as  the  ratio  of  the  orbital  frequency  to  the  inertial  frequency  multiplied 
by  the  azimuthal  wavenumber.  Based  on  observations  from  Hurricane  Gloria  (1985),  SM  showed 
that  the  square  of  the  local  Rossby  number  is  generally  less  than  unity  only  for  wavenumber  one 
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in  the  rapidly  rotating  region  of  an  intense  vortex.  A  goal  in  this  work  is  to  further  explore  the 
wavenumber  selection  mechanism  that  provides  the  basis  for  AB  theory. 

Tb  aid  in  understanding  the  fundamental  physical  processes  of  axisymmetrization,  previous 
works  using  quasi-geostrophic  shallow  water  models  and  nondivergent  barotropic  models  are  now 
reviewed.  Melander  et  al.  (1987)  studied  the  axisymmetrization  in  a  quasi-geostrophic  shallow  water 
model  using  apseudospectral  formulation.  The  model  was  initialized  using  finite-amplitude  elliptical 
vortices  with  modified  ‘tophat*  radial  profiles.  In  that  paper,  axisymmetrization  was  identified 
with  the  shedding  of  vorticity  filaments  outside  the  nearly  axisymmetric  vortex  core.  Rather  than 
decompose  the  results  into  radial  and  azimuthal  components,  the  vortex  evolution  was  diagnosed 
by  plotting  the  aspect  ratio  of  vorticity  isolines  as  a  function  of  time.  At  early  times  (t  less  than 
a  few  orbital  periods),  all  aspect  ratios  decreased  towards  unity  and  then  began  small-amplitude 
oscillations.  Melander  et  al.  identified  this  initial  transition  with  axisymmetrization. 

For  the  majority  of  cases  they  considered,  representation  of  dissipative  processes  at  small  scales 
was  parameterized  using  a  hyperviscosity  formulation  in  the  potential  vorticity  evolution  equation. 
This  could  conceivably  have  played  a  role  in  the  axisymmetrization  process.  However,  by  increasing 
the  value  of  the  hyperviscosity  used  in  their  model,  Melander  et  al.  convincingly  demonstrated  the 
invisdd  nature  of  the  axisymmetrization  process. 

In  a  complementary  approach,  McCalpin  (1987)  used  a  reduced-gravity  quasi-geostrophic  model 
to  study  axisymmetrization  in  Gulf  Stream  rings.  The  vortex  model  consisted  of  finite-amplitude 
azimuthal  mode  2  or  mode  3  perturbations  superposed  on  an  axisymmetric  Gaussian  basic  state. 
McCalpin  found  that  nearly  all  of  the  perturbation  energy  was  transferred  to  the  basic  state  on  time 
scales  on  the  order  of  an  orbital  period.  While  only  weakly  dependent  on  the  perturbation  strength, 
the  decay  time  scales  were  found  to  be  strongly  dependent  on  the  strength  of  the  mean  flow  and  the 
azimuthal  wavenumber.  In  particular,  the  decay  time  scale  for  wavenumber  three  was  found  to  be 
60%  faster  than  that  of  wavenumber  two. 

Like  Melander  et  al.,  McCalpin  parameterized  dissipative  processes  with  a  hyperviscosity  for¬ 
mulation.  However,  he  provided  only  a  limited  discussion  regarding  the  effects  of  diffusion  on  the 


2 


synunetrization  process.  McCalpin  calculated  a  diffusion  time  scale  that  appears  to  be  based  on 
L*/Ab ,  wliere  L  is  the  characteristic  horizontal  scale  and  Ab  is  the  hyperviscosity  coefficient.  For 
deformation  radius  scale  waves,  his  diffusion  time  scale  was  found  to  be  80  years.  However,  since 
the  differential  rotation  in  the  vortex  rapidly  reduces  the  radial  scale  of  the  asymmetries  (SM;  also 
Sutyrin,  1989),  one  must  not  neglect  this  effect  when  estimating  the  diffusion  time  scale. 

In  terms  of  McCalpin ’s  model  parameters,  an  estimate  of  the  diffusion  time  scale  that  incorpo¬ 
rates  the  differential  rotation  is 

U  ~  Afl[(47r2ro2n2+n2)/r2]2  (11) 

where  m  is  the  number  of  orbits  traversed  by  a  fluid  parcel  at  the  radius  of  maximum  winds  (r) 
and  n  is  the  azimuthal  wavenumber.  For  Ab  =  5  x  10®  m4s-1,  m  =  1,  n  =  2,  and  r  =  50  km,  tj  is 
found  to  be  approximately  13  hours.  In  calculating  id,  it  is  assumed  that  the  differential  rotation 
has  already  decreased  the  radial  scale  of  an  asymmetry  to  render  basic  state  quantities  effectively 
constant  in  a  first  approximation.  Based  on  the  results  of  Melander  et  al.,  this  occurs  in  roughly  a 
few  orbital  periods.  For  a  typical  run,  McCalpin  found  the  perturbation  decay  time  scale  to  be  1.5 
orbital  periods.  This  suggests  that  any  viscous  influence  was  likely  insignificant. 

Both  McCalpin  and  Melander  et  al.  simulated  asymmetric  vortex  evolution  using  fully  nonlinear 
numerical  models  that  generally  prevent  analytical  solutions.  While  their  results  represent  important 
and  meaningful  contributions  to  our  understanding  of  axisymmetrization,  basic  physical  processes 
may  be  masked  by  the  presence  of  nonlinear  dynamics,  diffusion,  and  the  /3-effect. 

Sutyrin  (1989)  developed  a  formal  solution  for  linear  disturbance  evolution  in  a  quasi-geostrophic 
shallow  water  model  on  an  /-plane.  In  this  model,  the  disturbance  potential  vorticity  was  conserved 
following  fluid  particles.  For  regions  where  the  basic  state  potential  vorticity  was  identically  zero, 
Sutyrin  showed  that  the  nonmodal  component  of  the  disturbance  potential  vorticity  became  oscil¬ 
latory  in  radius  with  the  oscillations  controlled  by  the  differential  rotation  of  the  fluid.  In  addition, 
the  radial  gradient  of  the  nonmodal  disturbance  potential  vorticity  was  shown  to  increase  linearly 
with  time.  Sutyrin  asserted  that  a  similar  nonmodal  disturbance  evolution  would  be  observed  in 
regions  with  a  continuous  and  monotonic  basic  state  potential  vorticity  profile. 
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The  linear  increase  with  time  of  the  nonmodal  disturbance  potential  vortirity  gradient  was 
fundamentally  due  to  the  reduction  of  the  radial  scale  of  the  disturbance.  Sutyrin  concluded,  then, 
that  the  symmetrization  process  was  due  to  the  differential  rotation  of  the  fluid  and  was  analogous 
to  perturbation  decay  in  rectilinear  simple  shear  flow  (Case,  1960). 

Like  Melander  et  al.,  Sutyrin  did  not  examine  the  possibility  of  an  invisdd  azimuthal  wavenum¬ 
ber  selection  mechanism  nor  did  he  investigate  the  energy  transfer  between  the  asymmetries  and  the 
circularly  symmetric  basic  state.  Identification  of  these  processes  is  complicated  in  Sutyrin’s  model 
since  the  solution  is  expressed  in  terms  of  an  infinite  function  space  operator. 

Carr  and  Williams  (1989,  hereafter  CW)  studied  axisymmetrization  using  an  inviscid  nondiver- 
gent  barotropic  vortex  model.  Their  evolutionary  model  consisted  of  small-amplitude  perturbations 
superposed  on  a  steady  axisymmetric  Rankine  flow.  CW  were  primarily  interested  in  the  asymmetry 
inducing  influences  of  (3  and  environmental  wind  shear. 

CW  asserted  that  the  damping  rate  of  perturbations  is  proportional  to  the  square  of  the  az¬ 
imuthal  wavenumber.  This  conclusion  was  based  on  results  derived  from  initial  conditions  in  which 
both  the  azimuthal  wavenumber  and  the  radial  structure  changed  simultaneously.  CW  asserted 
further  that  their  result  was  analogous  to  the  dependence  of  damping  on  the  zonal  wavenumber  as 
described  by  Case  for  plane  Couette  flow.  However,  a  careful  review  of  Case  did  not  reveal  any 
discussion  of  how  the  perturbation  damping  rate  and  the  zonal  wavenumber  are  related.  Since  this 
relationship  is  potentially  a  fundamental  aspect  of  an  inviscid  wavenumber  selection  mechanism  in 
sheared  flows,  both  Case  and  CW  must  be  revisited  to  clarify  this  issue. 

In  this  work,  the  axisymmetrization  process  will  be  studied  as  an  initial-value  problem  within 
the  framework  of  inviscid  nondivergent  vortirity  dynamics.  Since  the  curvature  vortirity  and  the 
curvilinear  coordinate  system  of  the  swirling  problem  may  complicate  understanding  of  the  distur¬ 
bance  dynamics,  the  analogous  rectilinear  simple  shear  problem  (Case;  also  Farrell,  1987)  will  be 
examined  first.  Identification  of  how  the  disturbance  evolution  depends  on  the  zonal  wavenumber 
and  the  meridional  structure  of  the  initial  condition  will  be  emphasized. 
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CHAPTER  2 


ALGEBRAIC  INSTABILITY  IN  UNBOUNDED  SIMPLE  SHEAR 


2.1  Introduction 

The  purpose  of  this  chapter  is  to  examine  the  evolution  of  small-amplitude  perturbations  in 
simple  shear  as  a  first  step  in  the  understanding  the  evolution  of  asymmetries  in  a  hurricane.  Of 
particular  interest  is  how  the  perturbation  evolution  depends  on  the  zonal  wavenumber  and  the 
meridional  structure  of  an  initial  vortidty  profile.  The  flow  is  assumed  unbounded  and  the  model 
employed  is  the  two-dimensional  invisdd  Euler  equations  on  an  /-plane.  Since  the  model  is  ul¬ 
timately  formulated  at  the  level  of  the  vortidty  equation,  the  Coriolis  parameter  does  not  affect 
interpretation  of  results,  but  is  retained  for  consistency  with  >  r  chapters. 

Linear  perturbation  theory  is  used  to  decompose  the  flow  into  a  meridionally  varying  basic  state 
and  small  perturbation  field.  Ultimately,  the  results  of  this  chapter  are  extended  to  hurricanes.  Since 
observations  show  that  the  asymmetric  winds  in  hurricanes  are  small  compared  to  the  symmetric 
tangential  wind  (SM),  the  use  of  linear  perturbation  theory  is  justified  on  observational  grounds. 

2.2  Rectilinear  Shear  Model 


The  zonal  and  meridional  momentum  equations  and  continuity  equation  in  Cartesian  coordi¬ 
nates  are,  respectively, 


du 

du 

du  , 

dt 

dv 

dv 

dv  , 

dt 

+”s 

+v»i+f — : 

du  dv 

_  i  „  _  ft 

dx  dy  ' 

dp 

dx' 

dp 


(2-1  a) 
(2-1  b) 
(2.1c) 


where  x  and  y  are  the  zonal  and  meridional  coordinates,  u  and  v  are  the  zonal  and  meridional  winds, 
/  is  the  constant  Coriolis  parameter,  and  p  is  the  pressure  divided  by  constant  density.  The  Coriolis 
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For  simple  shear  flow,  the  basic  state  wind  is  given  by  fi  =  Sy  where  5  is  a  constant  shear  (see 


figure  2.1).  The  basic  state  vorticity  is  then 


da 

c  =  _S. 

dy 


The  equations  of  motion  (2.6a)-(2.6c)  become 


dv!  dn'  .  ,  djf 

-  +  Sy—  +  f»=-- 

du1  dv'  n 
dx  +  dy 


while  the  vorticity  equation  (2.6d)  is 


(2.7) 


(2.8a) 

(2.86) 

(2.8c) 


(2.9) 


The  linearized  vorticity  equation  (2.9)  states  that  the  perturbation  vorticity  is  conserved  following 
the  basic  state  flow. 

For  notational  simplicity,  primes  are  now  dropped.  All  quantities  are  understood  to  be  per¬ 
turbations  unless  they  have  an  overbar.  Boundary  conditions  must  be  specified  before  (2.9)  can  be 
solved.  In  an  infinite  domain,  u  and  v  be  must  be  bounded  as  |x|  or  |y|  becomes  arbitrarily  large. 

To  satisfy  continuity  (2.8c),  a  perturbation  streamfonction  may  be  defined  such  that 


Equation  (2.3)  becomes 


C  =  vV 


(2.10) 


(2.11) 


where  the  Lapladan  operator  is  defined  by 


Equation  (2.9)  is  thus 


(2.12) 


(2.13) 
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There  are  several  methods  that  can  be  used  to  solve  (2.13).  The  normal  mode  and  initial-value 
approaches  are  considered  below. 

2.3  The  Normal  Mode  Approach 

In  the  normal  mode  approach,  separable  solutions  of  the  form 

#r,y,t)  =  *(  (2.14) 

are  assumed  where  k  is  a  real  zonal  wavenumber,  a  is  the  eigenfrequency,  and  $(y)  is  the  eigenfunc¬ 
tion.  Substituting  (2.14)  into  (2.13)  and  dividing  out  the  common  exponential  term  yields 

{cr  +  kSy)(j£s-k2'j*  =  0.  (2.15) 

If  the  factor  <r  +  kSy  is  naively  divided  out,  then  (2.15)  has  a  solution  of  the  form 

*  =  c1ekv  +  c2e-k*.  (2.16) 

The  boundary  conditions  for  (2.16)  require  that  9  be  bounded  as  y  -*  ±oo.  The  only  way  to  satisfy 
these  conditions  is  if  ci  =  ci  —  0.  Thus,  there  are  no  normal  modes  associated  with  a  +  kSy  ft  0 
and  the  naive  conclusion  is  that  there  are  only  trivial  solutions  to  the  problem.  Tb  see  why  this  is 
incorrect,  note  that  (2.15)  may  also  be  satisfied  if  a  +  kSy  =  0.  At  such  paints,  more  information  is 
needed  to  integrate  (2.15)  for  $. 

2.4  The  Initial- Value  Approach 

An  alternate  solution  strategy  follows  the  methodology  of  Case  (1960).  Since  the  basic  state 
flow  only  depends  on  y,  the  zonal  dependence  may  be  represented  by  a  Fourier  integral.  The  Fourier 
transform  pair  is  given  by 

OO 

/(*)=  J  f(x)e~ikxdx,  (2.17 a) 

— OO 

OO 

/(*)  =  £  J  f{k)e'kx  dk,  (2.175) 
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where  k  is  the  zonal  wavenumber.  Examination  of  the  initial- value  problem  associated  with  (2.13) 
is  of  interest.  This  problem  is  formulated  by  using  the  Laplace  transform  pair,  which  is  given  by 

OO 

f(9)(q)  =  f  Me-'1*,  (2.18a) 

o 

m=2 hifW{q)eqtdq'  {2186) 

where  q  is  the  Laplace  transform  parameter  (generally  complex)  and  C  denotes  a  line  parallel  to  the 
imaginary  axis  that  is  positioned  to  the  right  of  all  singularities  of  f^q\q). 

The  Fourier  transform  of  (2.13)  is 

+  =  0  (2.19) 

while  the  Laplace  transform  of  (2.19)  is 

<» + «*»>  -  **)  - **)  *-  <»■*» 

The  variable  4>0  is  the  meridional  structure  of  the  initial  streamfiinction  for  wavenumber  k  and  the 
right  side  of  (2.20)  is  the  initial  vortidty.  On  dividing  through  by  (q  +  ikSy),  equation  (2.20)  may 
be  written  as 

p-21> 

Noting  that  the  Laplace  transform  of  exp (—ikSyt)  is  l/(q  +  ikSy),  the  inverse  Laplace  transform  of 
(2.21)  is 

(J^-fc2)  rP  =  (2.22) 

The  boundary  conditions  require  that  ip  be  bounded  as  |y|  oo.  Recalling  the  a  parameter  from 
the  normal  mode  approach,  it  is  seen  that  <r  =  - kSy ,  or  a  +  kSy  =  0.  Unlike  the  separable  normal 
mode  solutions  which  maintain  a  constant  structure  in  time,  the  initial-value  problem  has  a  time 
dependent  meridional  structure. 

The  differential  equation  (2.22)  is  readily  solved  by  the  Green’s  function  method.  The  Green’s 
function  is  defined  as  the  solution  to 

( j? 2  ~  G(y»  Vo)  =  6(V  ~  Vo)  (2.23) 
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Equation  (2.31)  is  referred  to  as  the  continuous  spectrum  solution  to  equation  (2.13)  since  it  results 
from  integrating  over  all  advective  frequencies  -iki(y0)  spanned  by  the  flow. 

The  solution  is  completed  by  applying  the  inverse  Fourier  transform,  yielding 

OO  CO 

J  e'ks  J  G(y, y0)Co(fc> !/o)e_i*Svot dy0 dk.  (2.32) 

— OO  — OO 

In  this  chapter  and  the  next,  attention  is  restricted  to  initial  vortidty  perturbations  of  the  form 
f(y)eik’m  where  k,  is  a  pre-spedfied  wavenumber.  The  Fourier  transform  of  this  initial  condition  is 
f(y)6(km  -  k),  where  6(k,  —  k)  is  again  the  Dirac  delta  function.  When  an  initial  condition  of  this 
form  is  substituted  into  (2.32)  and  the  integration  is  performed,  the  delta  function  maps  every  k 
into  J k,.  For  simplidty,  the  delta  function  will  not  be  explidt  in  any  of  the  subsequent  manipulations 
and  k  will  be  used  to  refer  to  the  wavenumber  in  both  physical  and  spectral  space.  With  the  above 
convention,  the  solution  (2.32)  is  written  as 

OO 

Mx,V,t)=  J  G(y,Vo)Co(*, Vo)eifc(,~5vot)  dy0.  (2.33) 

—  OO 

In  (2.33),  the  subscript  it  is  a  reminder  that  the  delta  function  has  filtered  a  particular  wavenumber 
from  the  inverse  Fourier  transform  integral. 

Before  examining  the  specific  evolution  of  the  perturbation  streamfunction  for  various  initial 
vortidty  profiles,  recall  that  the  problem  is  invisdd.  Thus,  all  processes  are  time-reversible.  The 
convention  will  be  to  define  t  =  0  as  the  initial  time,  t  >  0  as  future  times,  and  t  <  0  as  past  times. 

2.5  The  Plane  Wave 


In  an  infinite  domain,  a  natural  choice  for  an  initial  condition  is  one  that  is  periodic.  Let 


■0o  =  ea°v 

where  le  is  the  initial  meridional  wavenumber.  The  initial  vortidty  is  then 

Co  =  — (fc3  +  ll)eUoV. 

Substituting  (2.35)  into  (2.33)  and  integrating  yields 

•A.  — _ k  +  Ip  cHkx+lU-kSt)v 1 

*2  +  (lo  -  kStf 


(2.34) 


(2.35) 


(2.36) 
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Taking  the  real  part  of  (2.36)  gives 


*k  =  p  «*(**  +  «o  -  kSt)v ]•  (2.37) 

The  quantity  l0  -  kSt  is  a  time  dependent  meridional  wavenumber,  l,  which  describes  the  ever- 
changing  amplitude  and  meridional  structure  of  the  perturbation.  The  meridional  wavenumber  and 
time  are  negatively  related  so  as  t  increases  l  decreases. 

As  an  aid  to  understanding  the  evolution  of  the  physical  fields,  consider  the  wave  vector 
rh  =  (k,l).  The  magnitude  of  the  wave  vector,  Jrrl),  is  the  total  wavenumber  of  the  perturbation. 
The  angle  the  wave  vector  makes  with  the  zonal  axis  (measured  in  a  counter  clockwise  sense  from 
the  positive  *  axis)  is  given  by  0  =  tan-1(l/fc).  The  wave  vector  is  perpendicular  to  the  perturbation 
contours,  so  0  ranges  from  -jr/2  to  ir/2. 

The  streamfunction  evolution  described  by  (2.37)  is  qualitatively  considered  for  fixed  k  and 
la.  The  maximum  streamfunction  amplitude  is  (k?  +  ll)/k 2  and  occurs  when  t  =  la/kS  ( l  =  0). 
Thus,  the  streamfunction  amplitude  grows  when  t  <  l0/kS  ( l  >  0)  and  it  decays  for  t  >  l0/kS 
(I  <  0).  When  the  streamfunction  amplitude  is  maximum,  the  streamfunction  field  is  parallel  to  the 
meridional  axis  and  the  number  of  waves  in  the  domain  is  given  by  k.  Since  t  and  l  are  negatively 
related,  increasing  t  results  in  decreasing  0,  so  the  streamfunction  field  rotates  clockwise.  However, 
increasing  (decreasing)  t  results  in  |nl|  decreasing  (increasing)  when  t  <  l0/kS  ( t  >  la/kS).  Thus, 
when  the  streamfunction  field  is  in  the  growth  (decay)  phase,  the  number  of  waves  in  the  domain 
decreases  (increases). 

Figure  2.2  shows  streamfunction  contours  for  k  =  l0  =  1  at  nondimensional  times  St  ~  -10.0, 
-1.0,  1.0,  and  10.0.  When  the  streamfunction  field  is  in  the  growth  phase,  the  contours  have  a 
negative  slope.  When  the  streamfunction  field  is  in  the  decay  phase,  the  contours  have  a  positive 
slope.  The  shear  for  this  problem  is  positive  and,  thus,  has  a  positive  slope  (figure  2.1).  Comparing 
the  slope  of  the  streamfunction  contours  with  the  slope  of  the  shear,  the  streamfunction  field  tilts 
against  the  shear  in  the  growth  phase  while  the  streamfunction  field  tilts  with  the  shear  in  the  decay 
phase.  When  the  streamfunction  field  is  upright,  its  amplitude  is  maximum. 

To  further  quantify  the  growth  and  decay,  an  expression  for  the  perturbation  kinetic  energy 
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(KE)  density  is  derived.  Substituting  (2.37)  into  (2.10)  gives  the  zonal  and  meridional  winds  for  the 
plane  wave  streamfunction.  These  are 


(l0-kSt)(k’  +  l*0) 
k~  k2  +  {la-kSt)2 

„  *(**  +  £) 

*  ~  k2  +  (/„  -  ksty 

The  KE  density  is  then  p( u2  +  v2)/2,  which  gives 

p  (k2+lo)2 

k  ~  2 [fc2  +  {lo  -  kSt)2] 


sin[fca:  +  (i0  -  kSt)y ], 
sinffci  +  (lo  -  kSt)y]. 


sin2[fcx  +  (la  -  kSt)y] 


(2.38a) 

(2.386) 


(2.39) 


where  p  has  been  set  to  one. 

Equation  (2.39)  gives  a  local  measure  of  the  KE  density.  A  more  useful  measure  which  isolates 
the  temporal  evolution  of  the  disturbances  is  an  integrated  KE.  For  the  plane  wave,  the  integration 
is  performed  over  one  wavelength  zonally  and  meridionally  since  the  integrated  KE  density  over  the 
entire  plane  is  infinite.  The  integration  yields 


,r>  (k2+ll? 

K  k'  ~  4(Jfc2  +  ( l0  -  kSt )2] 


(2.40) 


where  (£*)  is  the  KE  per  wave.  This  is  maximum  when  l0  -  kSt  =  0,  which  gives 

(k2  +  ll)* 


(Ek)mmx  — 


4  Jfc2 


(2.41) 


Thus,  the  normalized  KE  per  wave  (hereafter,  the  normalized  KE)  is 


(Ek) _ k2 

(Ek)  max  k2  +  ( l0  -  kSt) 2 ' 


(2.42) 


The  normalization  has  been  done  with  the  maximum  value  of  the  normalized  KE  rather  than  the 
initial  value  so  that  the  normalized  KE  ranges  from  zero  to  one.  from  (2.42),  the  instantaneous 
growth  (decay)  rate  is  found  to  be 


1  d{Ek)  __  2kS(l0  —  kSt) 

(Ek)  dt  ~  k2  +  (l0-kSt)2'  1  J 

Setting  the  time  derivative  of  (2.43)  to  zero  gives  the  maximum  growth  (decay)  rate  which  occurs 
at  St  =  =pl  +  l0/k.  This  corresponds  to  l  =  ±k,  0  =  ±ir/ 4,  and  a  normalized  KE  of  0.5  for  both 
times. 
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This  suggests  a  natural  evolution  time-scale  for  the  perturbation  energy  to  decay  to  half  its 
maximum  energy,  the  half-life  time.  For  the  plane  wave  perturbation,  the  half-life  time  is  1/5. 
Interestingly,  this  is  independent  of  k  and  l0 ■  Figures  2.3  and  2.4  show  the  normalized  KE  and  the 
growth  (decay)  rate  as  functions  of  time  for  k  =  1,  la  =  1;  k  =  1,  la  =  2;  and  k  =  2,  l„  =  1.  These 
figures  clearly  demonstrate  that  changing  k  and  l0  merely  shifts  the  energy  curve  along  the  time 
axis,  but  do  not  change  the  structure  of  the  energy  curve. 

The  half-life  time  for  sheared  disturbances  in  mid-latitude  cyclones,  hurricanes,  and  tornados 
may  be  estimated  as  follows.  Letting  U  and  L  denote  typical  horizontal  and  temporal  scales,  the 
half-life  time  is  ~  L/U.  In  moderate  mid-latitude  cyclones,  V  —  10  ms-1  and  L  ~  1000  km  giving 
a  horizonatal  shear  of  10~s  s-1  and  a  half-life  time  of  about  1  day.  In  the  eyewall  of  a  hurricane, 
[/  ~  50  ms-1  and  L  ~  50  km.  Thus,  the  horizontal  shear  is  10~3  s_1  and  the  half-life  time  is  about 
15  minutes.  Finally,  in  a  tornado,  U  ~  100  ms-1  and  L  ~  100  m  yielding  a  horizontal  shear  of  1  s-1 
and  a  half-life  time  of  1  second.  These  numbers  are  only  approximate,  but  they  serve  to  illustrate 
how  the  strength  of  the  shear  affects  the  evolution  time-scale. 

Perturbation  growth  and  decay  imply  an  energy  transfer  to  and  from  the  disturbance.  Under¬ 
standing  of  this  process  is  obtained  by  examining  the  perturbation  kinetic  energy  equation,  which 
is  formed  by  taking  the  dot  product  of  (u,v)  with  equations  (2.8a)  and  (2.8 b).  This  yields 


dE  _3E 
dt  +Udx 


+ 


d(up)  _  d(vp) 
dx  dy 


(2.44) 


Integrating  over  the  perturbation  domain  gives 


(2.45) 


When  the  perturbation  grows,  the  left  side  of  (2.45)  is  positive.  For  positive  shear,  the  momentum 
flux  term  is  negative  and  perturbation  zonal  wind  is  fluxed  equatorward.  The  effect,  then,  is  that 
the  shear  decreases  and  the  basic  state  provides  energy  for  the  perturbation’s  growth.  When  the 
perturbation  decays,  the  left  side  of  (2.45)  is  negative.  Here,  the  momentum  flux  is  positive  and 
perturbation  zonal  wind  is  fluxed  poleward.  The  effect  is  that  the  shear  increases  and  the  basic  state 
strengthens  at  the  expense  of  the  perturbation. 
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From  (2.45)  it  is  seen  that  the  perturbation  energy  tendency  is  a  nonlinear  process  since  it 
depends  on  uv.  Thus,  the  perturbation  energy  is  not  conserved  in  the  linear  problem  and  the  basic 
state  energy  is  viewed  as  an  infinite  reservoir.  It  can  be  shown  from  the  nonlinear  momentum 
equations,  though,  that  the  total  energy  in  this  system  is  conserved. 

For  the  remainder  of  this  work,  attention  is  restricted  to  disturbance  decay,  so  l0  =  0.  The 
equations  developed  in  this  section  become 


Co  =  -*2> 

(2.46) 

(k  =  —k2  cos[fc(x  -  Sty)], 

(2.47) 

tk  -  i  +  (5^2  cos[*(x  Sty)], 

(2.48) 

«*  =  i  +  sin(*(x  Sty)], 

(2.49) 

-k 

Vk  =  i  +  (St)2  sin[A:(x  Sty)], 

(2.50) 

Ek  -  i  +  (S<)2  sm2[fc(x  Sty)], 

(2.51) 

{■ Ek )  1 

{Ek)0  l  +  (St)2’ 

(2.52) 

1  d(E)  -2  S2t 

(Ek)  dt  1  +  (St)2' 

(2.53) 

The  upright  plane  wave  initial  condition  is  formally  equivalent  to  a  constant  initial  vorticity  profile. 
Thus,  perturbations  forced  by  either  constant  or  periodic  initial  conditions  in  simple  shear  flow  have 
global  energy  evolutions  that  are  independent  of  the  zonal  wavenumber  k. 

2.6  Gaussian  Initial  Condition 

The  plane  wave  considered  in  the  previous  section  provides  analytical  solutions  for  the  per¬ 
turbation  evolution  in  simple  shear.  However,  transient  atmospheric  forcings  are  more  realistically 
represented  as  isolated  disturbances.  To  simulate  a  short  wave  trough,  a  Gaussian  initial  vorticity 
profile  is  considered.  This  is  given  by 

=  (2.54) 

where  1/y/a  defines  the  characteristic  e-folding  length  scale  of  the  Gaussian.  Unlike  the  plane  wave 
initial  condition,  when  (2.54)  is  substituted  into  (2.33),  there  are  no  closed  form  solutions  available. 
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Using  the  trapezoidal  rule  to  evaluate  the  integral,  figures  2.5  and  2.6  plot  streamfunction  and 
vorticity  contours,  respectively,  for  the  initial  vorticity  profile  Co  =  exp(-y2).  The  rows  show  zonal 
wavenumbers  k  =  1,  2,  and  3,  respectively,  while  the  columns  denote  nondimensional  times  St  = 
0.0, 1.0,  2.0,  respectively.  There  are  several  differences  between  the  Gaussian  and  the  upright  plane 
wave  cases.  For  a  fixed  time,  the  streamfunction  amplitude  decreases  with  increasing  k.  In  addition, 
the  figure  suggests  that  that  higher  wavenumber  perturbations  decay  faster  than  lower  wavenumber 
perturbations.  In  the  upright  plane  wave,  the  streamfunction  amplitude  and  the  decay  rate  were 
independent  of  the  zonal  wavenumber.  Finally,  there  is  a  suggestion  that  for  large  time,  the  plane 
wave  solution  emerges. 

As  with  the  plane  wave  perturbation,  a  global  measure  of  the  KE  is  desired.  Here,  though,  the 
meridional  integration  is  from  —  oo  to  +oo.  Figure  2.7  shows  the  normalized  KE  as  a  function  of 
time  for  the  initial  vorticity  profile  Co  =  exp(-y2).  As  indicated  in  figure  2.5,  higher  wavenumber 
perturbations  decay  more  rapidly  than  lower  wavenumber  perturbations.  In  addition,  for  increasing  k 
or  t ,  the  normalized  KE  curves  asymptote  to  the  plane  wave  solution.  Interestingly,  the  wavenumber 
k  =  1  perturbation  decays  much  slower  than  all  higher  wavenumber  perturbations.  This  would  not 
be  surprising  in  a  viscous  fluid.  However,  in  the  above  results,  the  viscosity  is  identically  zero 
yet  higher  wavenumber  disturbances  still  decay  faster  than  smaller  wavenumber  disurbances.  This 
suggests  a  truly  inviscid  decay  mechanism  that  rapidly  damps  the  higher  wavenumbers  leaving  the 
lower  wavenumbers  behind. 

To  gain  further  understanding  of  the  zonal  wavenumber  dependence  in  the  evolution  of  merid¬ 
ionally  confined  initial  vorticity  disturbances,  it  is  instructive  to  examine  the  interaction  of  the  three 
terms  comprising  the  integrand  of  (2.33).  Figure  2.8  shows  the  Green’s  function,  G(y,y„);  the  phase 
function,  cos {kSty0);  their  product;  and  the  initial  vorticity  Co  —  exp(— y2)  as  functions  of  y0  for 
zonal  wavenumber  k  =  1.  Figure  2.9  is  the  same  plot  for  zonal  wavenumber  Jfe  =  3.  In  the  figures,  the 
Green’s  function  is  only  considered  at  y  =  0  and  the  -l/2k  factor  is  ignored  since  it  is  constant  in 
the  integration.  In  addition,  the  time  is  chosen  so  that  St  =  1.  As  k  increases,  the  Green’s  function 
narrows  while  the  phase  function  oscillates  more  rapidly.  The  effect  of  increasing  k,  then,  is  for  the 
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product  of  the  Green’s  function  and  the  phase  function  to  become  more  meridionally  confined.  For 
large  enough  fc,  the  meridional  scale  of  this  product  is  so  small  that  the  initial  vorticity  appears  ef¬ 
fectively  constant  in  the  integration.  This  is  reminiscent  of  the  upright  plane  wave  solution  in  which 
the  initial  vorticity  profile  is  constant  in  the  integration.  Thus,  while  the  meridional  scale  of  the 
initial  vorticity  generally  introduces  a  dependence  on  the  zonal  wavenumber  in  the  shear  dynamics, 
for  large  k  the  plane  wave  solution  should  emerge. 

For  increasing  k,  the  plane  wave  solution  may  be  recovered  from  a  meridionally  confined  initial 
condition.  This  can  be  demonstrated  analytically  as  follows.  The  streamfunction  amplitude  in 


Fourier  space  is 


OO 

V>(y)=  J  G(y,y0)Uvo)e-ikSt^  dVo. 


(2.55) 


As  k  becomes  large,  the  streamfunction  amplitude  can  be  approximated  by 


OO 

#(v)  »  Uv)  J  G(y,  y0)e~ikStVm  dya  1  <  k  <  oo. 


(2.56) 


This  integral  may  be  evaluated  exactly  giving  the  approximate  streamfunction  amplitude 


V,  -  _  --ikSty 

*  *2[l  +  b*)a] 


(2.57) 


Upon  applying  the  inverse  Fourier  transform  and  taking  the  real  part,  (2.57)  becomes 


,Pk~~  fc*[/+(St)2]  "  5<y)1’ 


(2.58) 


so  the  meridional  velocity  is 


(2.59) 


Now,  the  zonal  velocity  is 


„  -  <*C»(y)  coB[*(g  -  sty)]  _  Co  (y)st  ....  _ 
k~  dy  Jba(l  +  (5t)2]  k[l  +  (St)*}  ^  Sty^' 


(2.60) 


Assuming  k  large  enough  but  finite  such  that  the  first  term  on  the  right  side  of  (2.60)  is  small 
compared  to  the  second  term,  the  zonal  velocity  becomes 


(2.61) 
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The  KE  density  is  then  approximated  by 


Ek  = 


<2o(y) 


sin2[*(:r  -  Sty)] 


*2[1  +  (St)2] 

where  the  density,  p,  has  been  set  to  one,  while  the  integrated  KE  is  given  by 

OO 

* Ek)  =  2/fc3[l  4-  (St)2]  /  y)dV ‘ 

“OO 

Normalizing  (2.63)  by  its  initial  value 

OO 

(Ek)o  =  jp'  J  C l(y)dy, 

—  OO 

gives  the  normalized  KE  for  large  * 


(2.62) 


(2.63) 


(2.64) 


(Ek)  =  1 

(Ek)  o  1  +  (St)2‘ 


(2.65) 


This  is  the  normalized  KE  for  the  upright  plane  wave.  A  physical  interpretation  for  the  zonal 
wavenumber  dependence  is  presented  at  the  end  of  this  section. 

Figure  2.5  shows  the  plane  wave  solution  emerging  for  large  time.  As  time  increases,  the  phase 
function  oscillates  more  rapidly  while  the  Green’s  function  and  initial  condition  remain  unchanged. 
For  sufficiently  large  times,  the  half-wavelength  scale  of  the  phase  function  is  narrow  enough  so  that 
the  initial  condition  appears  locally  constant  which  is  analogous  to  the  upright  plane  wave.  Farther, 
double  integration  by  parts  of  (2.55)  subject  to  the  boundary  conditions  reveals  the  t~2  dependence 
of  the  streamfunction  amplitude  at  large  times. 

Figures  2.X0  and  2.11  plot  perturbation  streamfunction  and  vortidty  contours,  respectively,  for 
zonal  wavenumber  *  =  1.  The  rows  show  the  fields  for  initial  vortidty  profiles  Co  =  exp(-ay2) 
where  a  —  1,  2,  and  3,  respectively,  while  the  columns  denote  the  fields  at  nondimensional  times  St 
=  o.o,  l.o,  and  2.0,  respectively.  Figure  2.12  plots  the  corresponding  normalized  KE  as  a  function 
of  time  and  shows  that  decreasing  the  meridional  scale  of  the  initial  condition  results  in  slower 
energy  decay.  Since  the  energy  decay  is  related  to  the  half-wavelength  scale  of  the  phase  function, 
decreasing  the  scale  of  the  initial  condition  increases  the  time  required  for  a  given  amount  of  energy 
decay.  Thus,  small-scale  disturbances  asymptote  to  the  plane  wave  solution  more  slowly  than  large- 
scale  disturbances. 
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Physically,  disturbance  evolution  is  governed  by  the  conservation  of  disturbance  vorticity  so  the 
meridional  change  in  the  zonal  wind  across  a  disturbance  controls  how  rapidly  a  vorticity  disturbance 
is  elongated.  Because  the  area  of  a  fluid  element  does  not  change  in  nondivergent  flow,  an  elongating 
vorticity  patch  implies  an  attendant  decrease  in  the  disturbance  meridional  scale.  The  decrease  in 
the  meridional  scale  is  directly  related  to  the  decrease  in  disturbance  energy.  Disturbances  of  large 
meridional  extent  are  elongated  more  rapidly  than  those  of  small  meridional  extent  and,  thus, 
decay  more  rapidly.  Since  the  plane  wave  solution  spans  the  domain,  its  meridional  length  scale 
is  effectively  infinite.  Therefore,  the  plane  wave  always  decays  more  rapidly  than  any  meridionally 
localized  disturbance. 

In  a  similar  fashion,  the  dependence  on  the  zonal  wavenumber  may  be  explained.  For  an 
initially  upright  disturbance,  the  zonal  wavenumber  is  inversely  related  to  the  initial  zonal  extent  of 
a  disturbance.  Consider  two  disturbances  with  the  same  initial  energy  and  meridional  extent  where 
one  disturbances  has  twice  the  initial  zonal  extent  of  the  other.  In  a  given  amount  of  time,  the  zonal 
elongation  of  both  disturbances  is  the  same.  To  conserve  area,  the  change  in  the  meridional  scale 
of  a  disturbance  with  large  zonal  extent  is  less  than  that  of  a  disturbance  with  small  zonal  extent. 
Therefore,  a  lower  wavenumber  disturbance  decays  slower  than  a  higher  wavenumber  disturbance. 

2.7  Summary 

The  evolution  of  disturbances  in  unbounded  simple  shear  flow  has  been  examined  as  a  first 
step  in  understanding  the  evolution  of  asymmetries  in  hurricanes.  The  system  was  formulated  as  an 
initial  value  problem  and  two  types  of  initial  conditions  were  considered.  For  the  plane  wave  initial 
vorticity  profile,  the  disturbance  evolution  was  only  a  function  of  time  and  decayed  as  t-2.  Since 
transient  atmospheric  forcings  are  more  realistically  represented  by  isolated  disturbances,  Gaussian 
initial  vorticity  profiles  were  considered.  Unlike  the  plane  wave,  the  evolution  of  these  perturbations 
depended  on  both  the  meridional  scale  and  the  zonal  wavenumber  of  the  initial  condition.  The 
meridional  scale  determined  how  rapidly  the  disturbance  was  elongated.  Large-scale  disturbances 
elongated  and,  subsequently,  decayed  more  rapidly  than  small-scale  disturbances.  The  dependence 
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1 


of  the  evolution  on  the  zonal  wavenumber  illustrated  a  truly  in  viscid  mechanism  associated  with 
differential  advection  that  favors  the  decay  of  small-scale  disturbances  over  large-scale  disturbances. 
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Simple  Shear  Flow 


Fig.  2.1:  Simple  shear  flow.  Arrows  show  basic  state  zonal  wind  vectors,  dashed  line  shows  the 
shear  (5  =  1). 


<Ek>/<Ek> 
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Fig.  2.3:  Normalized  kinetic  energy  for  a  plane  wave  initial  vortidty  profile  with  k  =  1,  l0  =  1 
(dot);  k  =  1,  l0  =  2  (dash);  k  —  2,  l„  =  1  (dot  dash). 
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Fig.  2.4:  Growth  (decay)  rates  for  a  plane  wave  initial  vorticity  profile  with  k 
k  =  1,  l0  =  2  (dash);  k  =  2,  J0  =  1  (dot  dash). 
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Fig.  2.5  :  Perturbation  streamfunction  for  the  Gaussian  initial  vorticity  profile  Co  =  exp  (— y2) . 
Panels  (a)-(c)  show  contours  for  zonal  wavenumber  k  —  1  at  times  St  —  0.0, 1.0,  and  2.0,  respectively. 
Pamela  (d)-(f)  show  contours  for  k  =  2  and  panels  (g)-(i)  show  contours  for  k  =  3  at  the  same  times. 
The  contour  interval  is  3.76  x  10-2  m2s-1  in  all  panels.  Solid  contours  denote  tpk  >  0,  dotted 
contours  <  0. 
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Perturbation  Vorticity 
Infinite  Channel 


Fig.  2.6:  Perturbation  vortidty  for  the  Gaussian  initial  vorticity  profile  Co  =  exp  (-y2).  The 
panels  are  the  same  as  in  figure  2.5.  The  contour  interval  is  1.33  x  10-1  s_1  in  all  panels. 
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<Ek>/<Ek> 
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Normalized  Kinetic  Energy 
Infinite  Channel 


Fig.  2.7:  Normalized  kinetic  energy  for  the  Gaussian  initial  condition  Co  =  exp(-y2)  for  k  =  1 
(dot),  Jk  as  2  (dash),  and  k  =  3  (dot  dash).  The  solid  curve  is  the  normalized  kinetic  energy  for  the 
upright  plane  wave. 
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Integrand  of  Equation  (2.33) 
Infinite  Channel 
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Fig.  2.8:  The  integrand  of  equation  (2.33)  for  A;  =  1,  y  =  0,  and  St  —  1.0.  (a)  The  Green’s 
function,  (b)  The  phase  function,  (c)  The  product  of  (a)  and  (b).  (d)  The  initial  vortidty  profile 
Co  =  exp(-y2). 
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Fig.  2.10:  Perturbation  streamfunction  for  zonal  wavenumber  k  =  1.  Panels  (a)-(c)  show  contours 
for  the  initial  vortidty  profile  (0  =  exp  (~y2)  at  times  St  =  0.0,  1.0,  and  2.0,  respectively.  Panels 
(d)-(f)  show  contours  for  the  initial  vortidty  profile  Co  =  exp  {-2y2)  and  panels  (g)-(i)  show  contours 
for  the  initial  vortidty  profile  Co  —  exp  (-3y2)  at  the  same  times.  The  contour  interval  is  3.76  x  10-2 
m2s_1  in  all  panels.  Solid  contours  are  rpk  >  0,  dotted  contours  i>k  <  0. 
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Perturbation  Vorticity 
Infinite  Channel 


Fig.  2.11:  Perturbation  vorticity  for  zonal  wavenumber  k  =  1.  Panels  are  the  same  as  in  figure 
2.10.  Contour  interval  is  1.33  x  10-1  s-1  in  all  panels. 
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<Ek>/ <Ek> 


Normalized  Kinetic  Energy 
Infinite  Channel 


Fig.  2.12:  Normalized  perturbation  kinetic  energy  for  zonal  wavenumber  k  =  1  for  the  Gaussian 
initial  vortidty  profiles  Co  =  exp(-y2)  (dot),  Co  =  exp(-2y2)  (dash),  and  Co  =  exp(-3y2)  (dot 
dash).  The  solid  curve  is  the  normalized  kinetic  energy  for  the  upright  plane  wave. 
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CHAPTER  3 


ALGEBRAIC  DECAY  IN  BOUNDED  SIMPLE  SHEAR 


3.1  Introduction 

The  purpose  of  this  chapter  is  to  briefly  examine  how  meridional  boundaries  alter  the  evolution 
of  perturbations  in  simple  shear  flow.  Since  the  evolution  of  asymmetries  in  a  bounded  vortex  is 
the  topic  of  the  next  chapter,  investigating  the  effects  of  boundaries  in  the  simple  shear  problem  is 
a  natural  intermediate  step. 

3.2  Bounded  Rectilinear  Shear 


The  presence  of  boundaries  does  not  alter  the  governing  equations  used  in  chapter  2.  Neverthe¬ 
less,  meridional  boundaries  provide  a  convenient  length  scale  that  may  be  used  to  nondimensionalize 
the  equations.  Letting  (x,y)  =  L(x,  y),  (  =  S(,  and  t  =  t/S  where  L  is  the  boundary  separation 
and  the  tildes  denote  nondimensional  quantities,  the  linearized  vorticity  equation  (2.9)  becomes 


ac  M  n 

9t+Vdx  ~  • 


(3.1) 


Although  the  magnitude  of  the  shear  is  no  longer  explicit  in  (3.1),  the  evolutionary  time  is  scaled  by 
the  shear.  Upon  application  of  the  Fourier-Laplace  transforms  defined  in  chapter  2,  equation  (3.1) 
becomes 


(3.2) 


For  notational  simplicity,  the  tildes  have  been  dropped.  Unless  otherwise  noted,  all  quantities  are 
henceforth  nondimensional. 

Solutions  to  (3.2)  depend  on  the  choice  of  boundaries.  Here,  the  meridional  boundaries  are 
assumed  to  be  slippery  walls  at  y  =  a  and  y  =  b  where  a  <  b.  Thus,  v(a )  =  v(b)  =  0  to  satisfy 
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the  no  normal  flow  condition.  In  terms  of  &  Fourier  perturbation  streamfunction  amplitude,  the 
boundary  conditions  become  rj>{a)  =  i>(b)  =  0. 

The  Green’s  function  technique  is  used  to  solve  (3.2)  subject  to  the  boundary  conditions.  Im¬ 
posing  the  continuity  and  jump  conditions  from  the  previous  chapter  along  with  the  boundary 
conditions  yields 


G(v,y») 


_ 1 _ /  sinh[fc(y  -  a)]  sinh[fc(y0  -  6)],  o  <  y  <  ya 

Jksinh[fc(6  -  a)]  \  sinh[fc(y  -  &)]  sinh(fc(y0  -  a)j,  ya  <  y  <  b  ' 


(3.3) 


The  solution  to  (3.2)  is  then 


t> 

=  J  G(y,Va  Ko 


,—ikty0 


dya- 


(3.4) 


To  complete  the  derivation,  the  inverse  Fourier  transform  is  applied  to  (3.4)  giving 

b 

1>k=J  G(y,y0)C,eik^x~tVo)  dya, 

a 

where  the  subscript  A;  is  in  accord  with  the  convention  of  chapter  2. 


(3.5) 


3.3  Uniform  Initial  Condition 


The  initial  vortidty  given  by  (0  =  -k2,  which  corresponds  to  the  upright  plane  wave  of  chapter 
2,  is  considered  first.  Substituting  Co  into  (3.5)  yields 


i>k  = 


l  +  t2  l 


cos[A:(x  -  ty)] 


sinh[fc(j/  -  6)]  cos[fc(x  -  to)]  -  sinh[Jfc(y  -  o)]  cos[k(x  -  t6)]  )  tn  ^ 
sinh[jfe(6-o)]  J’  (  6) 


The  first  term  inside  the  braces  is  recognized  as  the  infinite  plane  wave  solution  while  the  second 
term  is  the  boundary  correction  that  ensures  i>k{o)  =  i>k{b)  =  0.  Figure  3.1  plots  streamfunction 
contours  for  initial  vortidty  C  =  ~k2  and  boundaries  at  o  =  -1,  5=1.  The  rows  show  k  =  1, 
2,  and  3,  respectively,  while  the  columns  denote  t  =  0.0,  2.0,  and  4.0,  respectively.  The  boundary 
correction  in  (3.6)  decreases  the  maximum  streamfunction  amplitude  relative  to  the  unbounded  case 
since  sinh[fc(y — 6)]  <  0  and  -  sinh[fc(i/-o)]  <  0  for  all  y,  while  sinh[fc(5-o)]  >  0.  As  k  increases  for  a 
given  boundary  separation,  the  magnitude  of  the  boundary  correction  decreases  and  the  maximum 
streamfunction  amplitude  increases  since  the  hyperbolic  sine  in  the  denominator  increases  faster 


34 


than  either  of  the  hyperbolic  sine  terms  in  the  numerator.  In  the  limit  k  — »  oo,  the  boundary 
correction  goes  to  zero  and  the  plane  wave  solution  is  recovered. 

The  evolution  of  perturbations  in  a  bounded  channel  forced  by  a  uniform  initial  vorticity  profile 
may  be  further  analyzed  by  examining  the  integrated  kinetic  energy  (KE)  associated  with  (3.6). 
This  is  given  by 

IF  \  -  *2  f(l  +  *3)  .  2{cos[fo(6  -  a)]  -  cosh[*(6  -  a)]} 

' 1  2(l  +  t»)2|  2  Jk(6-o)sinh[Jb(6-o)] 

_  cos(fct(b  -  o)]  co8h[fc(b  +  q)]{sinh(2A:6)  -  sinh(2fca)} 

2k(6  -  a)  sinh2[Jfc(6  -  a)] 

cos[fct(6  -  a)]  sinh[fc(6  +  a)]{co8h(2k6)  -  cosh(2ka)} 

'  +  2 k(b  -  a)  sinh2[i(6  -  o)j 

sinh[2fc(b  —  a)]  1 

2Jb(6  -  o)  sinh2[Jfc(b  -  o)]  J 

The  first  term  inside  the  braces  is  the  energy  of  the  infinite  plane  wave  while  the  other  terms  describe 
how  the  perturbation  energy  is  changed  by  the  presence  of  boundaries.  Since  the  first  term  decays  as 
t~2  while  the  remaining  terms  decay  as  t~4,  the  infinite  plane  wave  solution  emerges  for  increasing  t. 
Figure  3.2  shows  the  normalized  KE  as  a  function  of  time  for  initial  vorticity  Co  =  -fc2  for  boundaries 
at  o  =  -1,  6  =  1.  Unlike  the  infinite  plane  wave  case  where  the  energy  decay  was  independent  of 
k,  perturbations  decay  more  rapidly  with  increasing  wavenumber  for  a  constant  initial  condition  in 
a  bounded  channel.  Moreover,  the  k  =  1  perturbation  decays  significantly  slower  than  the  higher 
wavenumber  perturbations.  While  very  different  from  the  infinite  plane  wave,  disturbance  evolution 
for  a  uniform  initial  vorticity  profile  in  bounded  shear  is  qualitatively  similar  to  that  of  an  unbounded 
Gaussian  initial  condition. 

The  similarities  arise  because  each  problem  has  a  characteristic  length  scale.  In  the  unbounded 
problem,  the  initial  vorticity  profile  defined  the  length  scale.  For  a  constant  initial  condition  in 
a  bounded  channel,  the  boundary  separation  defines  the  length  scale.  From  (3.6),  it  is  seen  that 
changing  the  boundary  separation  for  a  fixed  k  is  equivalent  to  changing  k  for  a  fixed  boundary 
separation.  Thus,  as  6  -  a  increases,  the  boundary  correction  term  decreases  and  the  maximum 
streamfunction  amplitude  increases.  As  expected,  in  the  limit  of  infinite  boundary  separation,  the 
plane  wave  solution  is  recovered. 
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Similar  arguments  explain  the  effect  of  increasing  boundary  separation  on  the  perturbation 
evolution  as  described  by  (3.7).  As  b— a  increases,  the  denominators  of  the  boundary  terms  increases 
more  rapidly  than  the  corresponding  numerators  and  the  perturbation  decays  more  rapidly.  In  the 
limit  of  infinite  boundary  separation,  the  plane  wave  evolution  is  recovered.  Figures  3.1  and  3.2 
can  be  used  to  illustrate  the  described  behavior  if  the  plots  for  k  in  the  figures  are  reinterpreted  as 
plots  for  b  -  a.  Physically,  increasing  the  boundary  separation  increases  the  meridional  scale  of  the 
perturbation.  Thus,  the  perturbation  is  sheared  more  and  decays  more  rapidly.  This  is  analogous 
to  changing  the  scale  of  the  Gaussian  initial  condition  of  chapter  2. 

3.4  The  Gaussian 

The  evolution  of  perturbations  forced  by  Gaussian  initial  vorticity  profiles  in  a  bounded  channel 
is  now  considered.  Since  (3.5)  does  not  admit  closed  form  solutions,  the  trapezoidal  rule  is  used  to 
evaluate  the  integral. 

Figures  3.3  and  3.4  plot  streamfunction  and  vorticity  contours,  respectively,  for  initial  vorticity 
Co  =  exp (-y2)  and  zonal  wavenumber  Jfe  =  1.  The  rows  show  contours  for  boundary  separation  b-a 
—  0.2,  2.0,  and  20.0,  respectively,  while  the  columns  show  contours  for  time  t  =  0.0,  2.0,  and  4.0, 
respectively.  Note  that  the  bottom  row  is  virtually  identical  to  the  unbounded  case.  In  terms  of  the 
two  characteristic  length  scales,  the  rows  represent  b-a  <  l/y/a,  b  -  a  ~  1  /y/a,  and  6 — a  >  1  /y/a, 
respectively.  As  long  as  b  —  a  <  l /y/a,  the  maximum  streamfunction  amplitude  and  the  decay  rate 
will  be  decreased  relative  to  the  unbounded  case.  For  b-a>  1  /y/a,  the  unbounded  evolution  may 
be  recovered  provided  b  —  a  is  sufficiently  large  to  render  boundary  effects  in  the  Green’s  function 
negligible. 

Figure  3.5  shows  the  normalized  KE  as  a  function  of  time  for  initial  vorticity  Co  =  exp(— y2) 
and  zonal  wavenumber  k  =  1.  As  was  seen  for  the  constant  initial  condition  in  the  bounded  channel, 
increasing  the  boundary  separation  increases  the  decay  rate  since  the  perturbation  is  sheared  more. 
The  curve  corresponding  to  6— a  =  20.0  in  figure  3.5  is  the  decay  observed  in  the  unbounded  problem 
and  does  not  change  for  further  increases  in  the  boundary  separation.  At  this  point,  the  meridional 
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scale  of  the  perturbation  is  the  characteristic  length  scale.  For  this  curve  to  asymptote  to  the  plane 
wave  solution  at  small  times,  the  zonal  wavenumber  must  be  increased. 

3.5  Summary 

The  evolutionary  decay  of  perturbations  in  bounded  simple  shear  flow  has  been  examined  as 
an  initial-value  problem.  The  case  of  a  uniform  vortidty  profile  was  considered  first.  Relative  to 
the  unbounded  problem  of  chapter  2,  the  maximum  streamfunction  amplitude  was  smaller  and  the 
decay  rate  was  slower  because  the  presence  of  boundaries  introduced  a  characteristic  length  scale. 
The  infinite  plane  wave  is  not  a  localized  disturbance  and,  thus,  decays  faster  than  any  bounded 
disturbance.  In  the  limits  of  large  zonal  wavenumber,  boundary  separation,  and  time,  the  infinite 
plane  wave  solution  emerged. 

The  case  of  Gaussian  initial  vortidty  profiles  was  then  examined.  For  boundary  separations 
smaller  in  scale  than  the  meridional  scale  of  the  perturbation,  the  maximum  streamfunction  am¬ 
plitude  was  smaller  and  the  decay  was  slower  relative  to  the  unbounded  case.  The  scale  of  the 
boundary  separation  was  the  characteristic  length  scale.  When  the  boundary  separation  was  greater 
than  the  scale  of  the  perturbation,  the  unbounded  evolution  was  recovered  as  long  as  the  boundary 
separation  was  large  enough  to  make  the  boundary  effects  in  the  Green’s  function  negligible. 


37 


Fig.  3.1:  Perturbation  streamfunction  contours  for  initial  vortidty  £>  —  -&2  and  boundaries  at 
a  =  —1.0,  b  =  1.0.  Panels  (a)-(c)  show  contours  for  zonal  wavenumber  k  —  1  at  times  t  =  0.0,  2.0, 
and  4.0,  respectively.  Panels  (d)-(f)  show  contours  tor  k  =  2  and  panels  (g)-(i)  show  contours  tor 
Is  —  3  at  the  same  times.  The  contour  interval  is  6.21  x  10~3  in  all  panels.  Solid  lines  denote  V’t  >  0, 
dotted  lines  <  0. 
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Fig.  3.2:  Normalized  kinetic  energy  for  initial  vortidty  =  -k3  and  boundaries  at  a  =  —1.0, 
b  =  1.0.  The  curves  are  the  energy  decay  for  It  =  1  (dot),  k  =  2  (dash),  and  fc  =  3  (dot  dash).  The 
solid  curve  is  the  energy  decay  for  the  upright  plane  wave. 


Fig.  3.3:  Perturbation  streamfunction  contours  for  zonal  wavenumber  k  =  1  and  initial  vortidty 
Co  =  exp(— y2).  Panels  (a)-(c)  are  for  boundaries  at  a  =  -0.1,  b  =  0.1.  The  contour  interval 
is  3.43  x  10-4.  Panels  (d)-(f)  are  for  boundaries  at  a  =  —1.0,  b  =  1.0.  The  contour  interval  is 
2.11  x  10~3.  Panels  (g)-(i)  are  for  boundaries  at  a  =  -10.0,  b  =  10.0.  The  contour  interval  is 
3.76  x  10-2.  The  columns  denote  time  t  =  0.0,  2.0,  and  4.0,  respectively.  Solid  lines  are  V’t  >  0, 
dotted  lines  i>k  <  0. 
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Fig.  3.5:  Normalized  kinetic  energy  for  k  =  1  and  initial  vorticity  (0  =  exp(— y2).  The  curves 
are  the  energy  decay  for  boundaries  at  a  =  -0.1,  b  =  0.1  (dot);  a  =  —1.0,  b  =  1.0  (dash);  and 
a  =  —10.0,  6  =  10.0  (dot  dash) .  The  solid  curve  is  the  energy  decay  for  the  upright  plane  wave. 
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CHAPTER  4 


AXISYMMETRIZATION  IN  SWIRLING  FLOW 


4.1  Introduction 

The  purpose  of  this  chapter  is  to  study  the  axisymmetrization  process  in  a  vortex.  The  objective 
is  to  extend  the  results  of  the  previous  chapters  to  a  rapidly  rotating  fluid  with  radially  varying  shear. 
The  dependence  of  asymmetric  evolution  on  the  azimuthal  wavenumber  and  the  radial  structure  of 
an  initial  vortidty  profile  is  of  particular  interest.  As  in  previous  chapters,  the  model  employed  is 
the  inviscid  incompressible  Euler  equations  on  an  /-plane,  but  here  a  cylindrical  coordinate  system 
is  chosen  to  represent  the  problem.  Consistent  with  observations  (SM),  linear  perturbation  theory 
is  used  to  decompose  the  flow  into  a  circularly  symmetric  basic  state  of  swirl  and  small-amplitude 
asymmetries. 

In  formulating  the  asymmetric  vortex  model,  an  appropriate  coordinate  system  must  be  chosen. 
A  natural  coordinate  system  is  one  that  moves  with  the  center  of  the  vortex,  i.e.,  a  storm-relative 
coordinate.  Following  Willoughby  (1979),  SM  showed  that  in  storm-relative  coordinates,  vortex 
motion  on  an  /-plane  merely  modifies  the  definition  of  pressure  while  preserving  the  form  of  the 
momentum  equations.  Since  the  model  adopted  below  is  ultimately  formulated  at  the  level  of  the 
vortidty  equation,  the  storm  motion  is  no  longer  explidt  in  the  vortidty  dynamics  and  only  appears 
in  the  boundary  conditions  at  infinity.  In  a  Fourier  azimuthal  representation,  storm  motion  projects 
only  onto  the  azimuthal  wavenumber  one  component.  For  the  simple  case  of  a  resting  environment, 
the  assodated  winds  must  coindde  with  the  motion  at  infinity.  All  other  azimuthal  components 
must  vanish  at  infinity.  The  more  common  choice  of  coordinates  is  the  ground-based  system.  Here, 
vortex  motion  is  explidt  in  the  vortidty  dynamics  and,  for  a  quiescent  environment,  all  asymmetries 
must  vanish  at  infinity.  Unless  otherwise  stated,  a  ground-based  coordinate  system  will  be  adopted. 


43 


4.2  Swirling  Model 


The  radial  and  tangential  momentum  equations  and  continuity  equation  in  cylindrical  coordi¬ 
nates  are,  respectively, 


du  du  v  du  .  v 2  dp 

~dt  +Uar  +  rdX  ~*V~T  ~  ~!fr' 
dv  dv  v  dv  ,  uv  1  dp 

dt  dr  r  dX  r  r  dX 

1  d(ru)  1  dv  _ 
r  dr  rdX  ’ 


(4.1a) 


(4.1c) 


where  r  and  A  are  the  radius  and  azimuth,  u  and  v  are  the  radial  and  azimuthal  winds,  /  is  the 
constant  Coriolis  parameter,  and  p  is  the  pressure  divided  by  constant  density.  Explicit  dependence 
on  the  pressure  and  Coriolis  parameter  may  be  removed  by  talcing  the  horizontal  curl  of  (4.1a)  and 
(4.16).  The  result  is  the  vorticity  equation,  given  by 

Here  C  is  the  vertical  component  of  the  vorticity,  defined  by 

1  d(rv)  1  du 

^TaT-fax  (4'3) 

Equation  (4.3)  states  that  vorticity  is  conserved  following  the  fluid  particles. 

Since  intense  atmospheric  vortices  often  exhibit  weak  asymmetries,  linear  perturbation  theory 
is  applied  about  an  axisymmetric  basic  state  of  swirl.  Let  u  —  u’,  v  =  v (r)  +  v’,  and  p  =  p  +  jJ 
where  an  overbar  denotes  the  basic  state  and  a  prime  denotes  an  asymmetry.  Neglecting  products 
of  primed  quantities,  equations  (4.1a)-(4.1c)  become 


du ’  vdu’  .  ,  ,  v 2  2v  .  d  , 

«r + ;  ax  - {e  - ,v  ~  7  -  vv  “ >• 

dv'  vdv'  ,dv'  ,  ,  v  ,  1  3  , 

ar+;3A+"*+/“+;“=-;5A(j,+</)’ 


(4.4a) 


d(ru')  dv' 

-ir+a x=a  (4-4c) 

Noting  from  (4.3)  that  the  basic  state  vorticity  is  only  a  function  of  r,  the  vorticity  equation  (4.2)  is 


dC  vd('  (dC  „ 

+  ~H\  +  u  ~T  ~ 
dt  r  d  X  dr 


The  quantities  fj  and  £  represent  the  absolute  vortidty  and  the  effective  Coriolis  parameter  of  the 
symmetric  vortex.  For  notational  simplicity,  the  primes  will  now  be  dropped.  All  quantities  will  be 
perturbations  unless  they  have  an  overbar.  Attention  is  now  restricted  to  (4.8d). 

Equation  (4.8d)  is  nondimensionalized  by  defining  r  =  Rmf,  (u,v)  =  Vm(u,v),  t  =  Rm/Vmt, 
(C»*l)  =  and  A  =  A  where  Vm  is  the  maximum  azimuthal  wind,  Rm  is  the  radius  of 

maximum  winds,  and  tildes  denote  dimensionless  quantities.  Substituting  these  definitions  back  into 
equation  (4.8d),  the  factor  V%/Rm  is  found  common  to  every  term.  On  dividing  out  this  term,  the 
nondimensional  vortidty  equation  is  formally  identical  to  equation  (4.8 d).  Dropping  the  tildes  for 
notational  simplidty,  the  nondimensional  linear  vortidty  equation  is 


dC  t;  dQ  dfj 

^  +  rTx+UTr 


=  0 


(4.9) 
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Unlike  the  rectilinear  problem  in  simple  shear  where  there  was  no  gradient  of  the  basic  state  absolute 
vortidty,  the  perturbation  vortidty  in  a  swirling  flow  is  generally  not  conserved  following  the  basic 
state  flow. 

For  the  rectilinear  problem  in  simple  shear,  the  meridional  gradient  of  the  basic  state  zonal  wind 
controlled  the  perturbation  decay.  The  vortidty  tendency  was  directly  related  to  the  advection  of 
the  vortidty  by  the  basic  state  zonal  wind.  For  the  swirling  problem,  the  vortidty  tendency  is 
proportional  to  the  advection  of  the  vortidty  by  the  symmetric  vortex.  The  basic  state  angular 
velodty  is  playing  an  analogous  role  in  the  swirling  environment  that  the  basic  state  zonal  wind 
played  in  the  rectilinear  environment.  It  is  therefore  anticipated  that  the  radial  gradient  of  the  basic 
state  angular  velodty  will  control  the  decay  of  asymmetries  in  the  swirling  vortex. 

Before  equation  (4.9)  can  be  solved,  boundary  conditions  must  be  spedfied.  Since  observations 
show  that  the  rotation  of  a  hurricane  inside  the  radius  of  maximum  winds  is  in  approximate  solid 
body  rotation,  the  angular  velodty  gradient  is  small  in  that  region  and  little  symmetrization  is 
expected  there.  To  isolate  the  basic  decay  process,  a  boundary  is  placed  at  the  radius  of  maximum 
winds.  While  this  approach  is  justified  on  observational  grounds,  it  also  has  the  mathematical 
advantage  of  avoiding  .he  geometric  singularity  assodated  with  the  vortex  center.  Since  observations 
also  show  that  asymmetries  become  comparable  to  the  symmetric  azimuthal  wind  several  hundred 
kilometers  from  the  storm  center,  the  linear  assumption  is  adequately  satisfied  for  an  outer  boundary 
placed  so  that  v  >  («',«')•  In  nondimensional  coordinates,  r  —  a  denotes  the  boundary  at  the  radius 
of  maximum  winds  while  r  =  6  gives  the  location  of  the  environmental  boundary.  The  boundaries 
are  assumed  to  be  slippery  walls  at  which  the  normal  velocity  vanishes.  Thus,  u(a)  =  u(b)  =  0. 

To  ensure  that  solutions  to  the  vorticity  equation  (4.9)  satisfy  continuity  (4.8c),  define  a  per¬ 
turbation  streamfunction  such  that 


Equation  (4.3)  becomes 


1  dip 

^  A  %  J 

r  a  X 


C  =  vv 


(4.10) 


(4.11) 
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where  the  Laplacian  operator  is  given  by 

v’=f— (V-V1— L 

\rdr\dtj  rJ0A>J 

The  vorticity  equation  (4.9)  is  then 

(d  vd\\\d(d\  1  a2 1  ,  i  df}di>  n 

(at  +  raAj  [rdr  (r0r  J  +  r2  dA2J  ^  r  dr  dX  ~  '  *412) 

Equation  (4.12)  describes  the  evolution  of  the  system  subject  to  the  constraint  that  the  az¬ 
imuthal  gradient  of  the  streamfunction  vanish  on  the  boundaries.  Further  transformation  of  this 
equation  is  desired  to  facilitate  finding  its  solutions.  Because  the  azimuthal  variable  is  periodic 
by  definition,  the  azimuthal  dependence  may  be  represented  by  a  Fourier  series  where  the  Fourier 
transform  pair  is 

2t 

/(«)  =  J  f(X)e~inX  dX;  n  =  ±1,  ±2,  ±3, . . .  (4.13a) 

o 

/( A)=  J  /(»)e*nA;  »*0  (4.136) 

n=— oo 

and  n  is  the  azimuthal  wavenumber.  The  wavenumber  zero  contribution  is  excluded  from  the  above 
definitions  because  it  represents  the  symmetric  component  which  does  not  change  in  the  linear 
problem.  The  Fourier  transform  will  hereafter  be  applied  to  asymmetric  quantities.  As  in  chapter  2, 
the  initial-value  problem  associated  with  (4.12)  is  considered.  This  is  done  by  using  the  Laplace 
transform  pair,  given  by 


<30 

/(,)(«)  =  Jf(t)e-«dt 
o 

m  =  - h£f{9){q)e,tdq 


(4.14a) 

(4.146) 


where  q  and  C  are  the  same  as  defined  in  chapter  2. 

The  Fourier  transform  of  (4.12)  is 

/  d  V\  fl  d  (  d  \  n2l  ;  indfj  7 

while  the  Laplace  transform  of  (4.15)  is 


(4.15) 


(’  +  “')  [ll  (rl)  '  £]  *f"  ~  7**“  =  [r  J?  (r|)  -  S]  <416> 
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with  i>0  the  radial  structure  of  the  initial  streamfunction  for  wavenumber  n.  The  right  side  of 
equation  (4.16)  is  the  initial  vortidty  for  wavenumber  n.  Upon  dividing  by  q  +  inv/r,  this  equation 
can  be  rewritten 

p*frM  ga1  .;•(«>  .  =  4  (4.17) 

[rdr  \  dr)  r2J  r  dr  lq  + inv/r)  (q  +  in€/r) 

Unlike  the  rectilinear  problem  with  zero  vorticity  gradient,  the  inverse  Laplace  transform  of 
(4.17)  does  not  simply  result  in  a  differential  equation  for  V>  (cf.  (2.22)).  Rather,  the  direct  inversion 
of  the  second  term  on  the  left  side  of  (4.17)  requires  the  convolution  theorem  and  yields  tp  as  the 
argument  of  a  time  integral.  An  alternative  solution  method  for  the  analogous  rectilinear  problem 
is  described  by  Case  (1960)  which  requires  evaluation  of  a  contour  integral  in  the  complex  plane 
involving  a  Green’s  function  for  (4.17)  multiplied  by  the  right  side.  Ultimately,  numerical  schemes 
are  needed  to  obtain  explicit  solutions.  Such  procedures,  and  thus  the  general  solution  to  (4.17),  are 
beyond  the  scope  of  this  thesis. 

Here,  attention  is  restricted  to  a  specific  class  of  vortices  that  admit  simple  solutions.  If  the 
second  term  on  the  left  side  of  the  equation  is  zero,  the  solution  to  (4.17)  is  readily  constructed.  For 
the  geometry  considered  in  this  chapter,  the  Rankine  profile  (see  figure  4.1)  given  by 

v  =  o/r,  r  >  a  (4.18) 


satisifes  this  requirement  since  fj(r)  =  f.  Equation  (4.18)  then  becomes 


[i£(4)-£ 

[r  or  \  or  J  r2 


xpW  = 


(4.19) 


q  +  ian/r 2 ' 

Noting  that  the  Laplace  transform  of  exp(-ian</r2)  is  l/(q+ian/r2),  the  inverse  Laplace  transform 
of  (4.19)  yields 

[;£('!) <420> 


The  Green’s  function  method  will  now  be  used  to  solve  (4.20).  The  equation  is  first  put  in  standard 
form  by  multiplying  by  r,  giving 


=  rC0e-,ant/rJ. 


(4.21) 
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It  is  easily  verified  that  G(r,p)  =  G(p,r).  Thus,  in  an  analogous  manner  to  chapter  2,  the  solution 
to  (4.21)  is 

b 

i>(r)  =  J  G(r,p)Up)e-tant/p'  pdp.  (4.27) 

a 

In  the  previous  chapters,  only  initial  vorticity  profiles  of  the  form  Co  =  /(y)6(fc.  —  k)  were  con¬ 
sidered.  This  allowed  examination  of  the  perturbation  evolution  on  a  wavenumber  by  wavenumber 
basis.  A  similar  strategy  is  employed  here.  The  initial  vorticity  will  be  restricted  to  functions  of 
the  form  Co  =  /(r)tf(n.  -  n)  where  6  is  the  Kronecker  delta.  The  effect  on  the  integration  will  be 
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the  same  as  in  the  previous  chapters.  The  6  maps  n  into  n.  when  the  inverse  Fourier  transform  is 
applied.  For  notational  simplicity,  the  variable  n  is  used  to  refer  to  the  wavenumber  in  both  physical 
and  Fourier  space.  With  this  convention,  the  inverse  Fourier  transform  of  equation  (4.27)  is 

k 

1>n  =  emA  J  G(r,  pdp.  (4.28) 

a 

Since  exact  solutions  to  (4.28)  cannot  generally  be  found,  the  trapezoidal  rule  is  used  as  the  solution 
method  throughout  the  remainder  of  the  chapter.  Note  also  that  the  initial  condition  in  (4.28)  is 
expressed  in  terms  of  the  initial  vortidty.  Hereafter,  the  terms  initial  vortidty  and  initial  condition 
will  be  used  interchangeably.  Finally,  the  locations  of  the  boundaries  are  set  at  a  =  1  and  b  =  10. 

4.3  Initial  Conditions 

In  order  to  gain  insight  into  the  asymmetric  dynamics  that  are  applicable  to  hurricanes  and 
other  rapidly  rotating  vortices,  initial  vortidty  profiles  are  needed  that  represent  physical  processes 
relevant  in  such  vortices.  CW  motivate  four  classes  of  initial  conditions.  Their  analysis  is  outlined 
below.  Consider  two-dimensional  incompressible  invisdd  vortidty  dynamics  on  a /9-plane.  In  ground- 
based  coordinates,  the  governing  equation  is 

|  +  ^-V((  +  /)  =  0  (4.29) 

where  £  is  the  vertical  component  of  the  vortidty,  V*  is  the  horizontal  wind  vector,  and  /  is  the 
Coriolis  parameter  given  by 

f-fo  +  0y-  (4.30) 

Here  0  is  the  linearized  latitudinal  derivative  of  the  Coriolis  parameter  at  a  reference  latitude  and 
fo  is  the  reference  value  of  /.  Equation  (4.29)  may  be  transformed  into  a  coordinate  system  moving 
with  the  storm  center  through 

V  =  V',  (4.31) 

|  V',  (4.32) 

where  C(t)  is  the  storm  motion  vector  and  the  primes  denote  operations  with  respect  to  coordinates 
in  the  moving  frame. 
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The  velocity  is  now  partitioned  into  three  parts 

V  *  +  £  +  Va  (4.33) 

subject  to  the  following  definitions:  (a)  V,  is  a  known  symmetric  flow  that  is  steady  in  the  moving 
reference  frame,  (b)  V',  is  a  known  zonal  environmental  flow  that  depends  only  on  latitude  and  is 
steady  in  the  stationary  reference  frame,  and  (c)  Va  is  an  asymmetric  flow  that  represents  unknown 
perturbations  to  the  symmetric  flow.  Using  these  definitions,  the  vorticity  equation  in  the  moving 
reference  frame  is 

+  (?.  +  V.  -  d)  •  v'c,  +  (K  -  d)  •  V'c,  +  ?a  •  V'(C,  +  c«  +  f) 

=  -^V'(C,  +  /')-VeV'C..  (4.34) 

The  terms  have  been  grouped  so  that  external  processes  which  act  to  generate  Co  appear  as  forcing 
terms  on  the  right  side  of  (4.34). 

The  initial  vorticity  perturbations  are  obtained  by  examining  the  forcing  terms  in  component 
form.  CW  define  the  vectors  as  follows: 


=  Uj(r)A, 

(4.35a) 

=  u,(y)i, 

(4.356) 

<?  =  C[cos(A  -  a)r  —  sin(A  -  a) A], 

(4.35c) 

where  t,  r,  and  A  are  the  zonal,  radial,  and  azimuthal  unit  vectors,  respectively.  The  angle  a  is 
the  direction  of  the  storm’s  motion  following  the  convention  of  Willoughby  (1988).  CW  ote  that 
asymmetries  can  be  forced  by  the  term  -C  •  V'C,  which  is  the  advection  of  the  symmetric  vorticity 
by  the  storm’s  motion.  Grouping  this  term  with  the  other  forcing  terms,  the  vorticity  tendency  in 
Cartesian  coordinates  due  to  these  forcings  is 

— 1  «  -v,  (~fy+P)  cosA  ~  (««  -c^^-cosA  +  Cy^sinA  (4.36) 

where  c*  =  Ceos  a  and  Cy  =  Csin  a.  To  assess  the  radial  structure  of  the  forced  response,  consider 
the  truncated  Taylor  series  expansion  of  u«  about  the  position  of  storm’s  center 

«e(y)  =  «e(0)  +  Sey  -  1 0eV 2  (4.37) 
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where 


s*  ~  1  , 

ay  t  v— o 


&  =  - 


dy 3  I »=o 


(4.38a) 

(4.386) 


The  scaling  for  v,  and  <K»/dr  may  be  identified  by  examining  the  modified  Rankine  profile 

x 


The  radial  derivative  of  C,  is  then 


_  (V2  .  ,  VmRm 

dr  ~  '  rx+2  ' 


(4.39) 


(4.40) 


The  Rankine  profile  used  in  this  chapter  is  recovered  by  setting  X  =  1  in  (4.39).  In  order  to  have 
a  nontrivial  scaling  on  (4.40),  however,  X  cannot  be  exactly  one.  It  is  assumed  that  for  X  very 
close  to  one,  the  vortidty  gradient  will  be  sufficiently  small  that  the  perturbation  evolution  will  be 
essentially  unaffected.  Thus,  equation  (4.40)  becomes 


0C*  /  vJ  1  \  VmRm 

dr  ~  }  r3  ' 


(4.41) 


Examining  the  first  term  on  the  right  side  of  (4.36),  substituting  (4.386)  and  (4.39)  gives 


?k'K-VnSz.u)+fl')ax\. 


(4.42) 


This  represents  the  generation  of  wavenumber  one  asymmetry  due  to  the  advection  of  absolute 
environmental  vortidty  (the  0  terms)  by  the  symmetric  vortex.  CW  refer  to  this  as  the  ’^-induced 
asymmetry’.  An  equivalent  nondimensional  intial  condition  is  then 


Co  =  n  =  l. 


(4.43) 


Examining  the  second  term  on  the  right  side  of  (4.36),  substituting  the  first  two  terms  of  (4.37) 
and  (4.41)  gives 

d<U 


9t 


oc  (1  -  X^VnRn  cos  A  +  sin  2A 


(4.44) 
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where  the  identities  y  =  r  sin  A  and  2  sin  A  cos  A  =  sin2A  have  been  used.  This  represents  the 
generation  of  a  wavenumber  one  asymmetry  due  to  the  vortex  motion  relative  to  the  environmental 
wind  and  a  wavenumber  two  asymmetry  due  to  the  linear  shearing  of  the  symmetric  vortex  by 
the  relative  environmental  wind.  CW  refer  to  these  as  the  ‘motion-induced’  and  ‘shear-induced’ 
asymmetries,  respectively.  Equivalent  nondimensional  initial  conditions  are  then 

C.«£.  * =  1;  (4.45) 

Co  =  n  =  2.  (4.46) 

The  last  initial  condition  used  by  CW  was  motivated  by  a  desire  to  have  exact  solutions  to 
(4.28).  By  making  the  substitution  p  =  1/p3,  it  can  be  shown  that  an  initial  vortidty  profile  of  the 
form 

Co  =  l/rn+4  (4.47) 

admit  analytical  solutions  to  the  integration  (cf.  equations  (30)  and  (31)  in  CW).  These  initial 
vortidty  perturbations  are  radially  confined  to  the  inner  boundary  and  may  be  physically  interpreted 
as  asymmetries  forced  by  convection  in  the  eyewall  of  tropical  cydones.  CW  refer  to  these  as 
‘convection-induced’  asymmetries. 

4.4  Model  Results 

Since  one  of  the  objectives  of  this  thesis  was  to  address  some  of  the  unresolved  issues  of  the 
CW  paper,  a  necessary  first  step  was  to  be  able  to  reproduce  their  results.  Although  their  stream- 
function  fields  were  easily  reconstructed,  the  integrated  kinetic  energy  plots  associated  with  the  four 
initial  conditions  developed  in  the  previous  section  were  not  reproduced  exactly.  Figure  4.2  is  a 
reproduction  of  figure  4a  from  CW  while  figure  4.3  shows  the  results  from  this  work.  Although  the 
trends  are  the  same,  the  magnitudes  and  shapes  of  the  curves  are  different.  Specifically,  CW  show 
a  more  rapid  energy  decay  for  ‘convection-induced’  initial  conditions  and  a  more  abrupt  transition 
to  the  asymptotic  behavior.  The  conclusions  CW  draw  from  their  figure,  however,  are  essentially 
qualitative  in  nature  and  can  still  be  obtained  from  figure  4.3. 
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CW  provide  little  information  about  the  numerical  model  used  to  generate  their  results.  It  is 
likely,  then,  that  the  conditions  of  their  experiments  have  no.  oeen  exactly  duplicated.  However,  we 
are  confident  that  the  results  presented  in  this  thesis  are  robust.  The  numerical  results  here  have 
been  verified  to  satisfy  the  momentum  and  mass  conservation  equations  to  within  the  discretization 
error  associated  with  the  trapezoidal  rule.  Moreover,  as  the  model  resolution  was  increased,  the 
results  were  quadratically  convergent  (see  Appendix  B). 

In  their  paper,  CW  restrict  attention  to  the  initial  vorticity  profiles  developed  in  the  previous 
section.  While  physically  motivated,  they  change  both  the  radial  structure  and  the  azimuthal 
wavenumber  simultaneously.  Using  the  CW  initial  conditions  as  a  reference,  initial  vorticity  profiles 
where  only  one  parameter  changes  at  a  time  are  considered. 

Figure  4.4  shows  the  streamfunction  fields  for  an  upright  initial  vorticity  distribution  Co  =  1/r3 
and  figure  4.5  shows  the  corresponding  vorticity  fields.  The  rows  represent  wavenumbers  n  =  1, 
2,  and  3,  respectively,  while  the  columns  designate  times  t  —  0.0,  3.6,  and  7.2,  respectively.  These 
figures  show  that,  with  respect  to  increasing  the  time  or  the  wavenumber,  the  rectilinear  and  swirling 
problems  behave  similarly.  This  is  not  surprising  as  the  Green’s  function  and  the  phase  function  are 
formally  identical  between  the  rectilinear  and  swirling  cases. 

While  the  mathematical  solutions  are  qualitatively  similar  between  the  rectilinear  and  swirling 
cases,  there  are  physical  differences  in  the  disturbance  evolution  between  the  two  systems.  Figure 
4.6  shows  the  energy  decay  in  the  rectilinear  problem  for  the  initial  condition  Co  =  -l/(y  —  ll)3. 
Figure  4.7  shows  the  energy  decay  in  the  swirling  problem  for  Co  =  1/r3.  These  plots  show  that,  on 
a  wavenumber  by  wavenumber  basis,  the  energy  decays  more  slowly  in  the  swirling  problem  than  in 
the  rectilinear  case.  In  addition,  the  figures  suggest  that  the  difference  between  the  wavenumber  one 
disturbance  and  the  higher  wavenumber  disturbances  is  more  pronounced  in  the  swirling  problem. 
These  differences  are  explained  by  examining  how  the  shear  acts  in  each  system. 

In  the  rectilinear  problem  considered  in  chapter  2,  the  shear  is  constant  an  defines  limiting 
values  for  the  energy  decay  of  the  perturbation.  In  particular,  the  limiting  energy  decay  is  given  by 
1/[1 4-  (St)3].  Thus,  the  limiting  time  for  a  perturbation  to  decay  to  half  its  initial  energy  is  t  =  1  /S, 


the  half-life  time.  In  the  swirling  problem,  however,  the  shear  varies  radially  and  is  only  a  local 
measure  of  the  flow.  In  order  to  investigate  the  limiting  integrated  energy  decay  of  asymmetries,  a 
global  measure  of  the  shear  is  required. 

In  simple  shear  and  in  the  bounded  Rankine  vortex,  the  physical  process  governing  disturbance 
evolution  is  the  conservation  of  vorticity  following  the  basic  state  flow.  In  the  rectilinear  problem  of 
chapter  2,  the  perturbation  vorticity  is  advected  by  fi,  the  basic  state  zonal  wind.  The  corresponding 
shear  is  thus  du/dy.  Since  the  shear  is  constant,  the  meridional  scale  of  the  perturbation  qualitatively 
describes  the  decay  rate  for  a  constant  zonal  wavenumber.  The  change  in  the  zonal  wind  across  a 
large-scale  perturbation  is  greater  than  the  corresponding  change  across  a  small-scale  perturbation. 
Thus,  the  large-scale  disturbances  decay  more  rapidly. 

In  the  swirling  problem,  the  disturbance  vorticity  is  advected  by  v/r,  the  basic  state  angular 
velocity.  The  corresponding  local  shear  is  rdfl/dr  (SM)  where  Cl  =  C/r.  In  a  Rankine  vortex.  Cl  is 
equal  to  1/r2  and  the  local  shear  is  equal  to  — 2/r2.  Since  the  angular  velocity  does  not  describe 
solid  body  rotation,  both  the  radial  location  and  the  radial  structure  of  the  asymmetry  are  needed 
to  describe  the  decay  rate.  A  plausible  candidate  for  an  effective  shear  that  satisfies  these  conditions 
is 


Ho  dr 

a 

The  analogous  expression  in  the  rectilinear  problem  reduces  to  the  simple  shear  definition  if  du/dy 
is  constant.  For  initial  conditions  of  the  form  l/r“,  the  effective  shear  is  found  to  be 


5eff  =  - 


2(q  -  1)  f  oa+1  -  ba+l 


Stfl  = 


abfiia  + 1) 
a2- 6s 
062  ln(6/a)  ’ 


+1  -  6a+l' 

-1  ’ 


a  ±  1; 


a  =  1. 


(4.49a) 

(4.496) 


Appealing  to  the  rectilinear  case,  the  effective  shear  should  define  limiting  values  for  the  decay  rate 
of  asymmetries.  In  particular,  the  limiting  energy  decay  should  be  given  by  1/[1  +  (Seffi)2].  Thus, 
the  limiting  asymmetry  half-life  time  should  be  t  =  —  1  /S,tt-  The  limiting  energy  decay  curve  is  the 
solid  line  in  figure  4.7. 


55 


Figure  4.8  shows  the  integrated  energy  as  a  function  of  time  for  azimuthal  wavenumber  n  =  1 
and  vortirity  profiles  £,  =  l/ra  where  a  =  3,  4,  5,  and  6.  Figure  4.8  confirms  the  expected  behavior 
for  the  1/r3  and  1/r4  asymmetries.  However,  the  decay  rates  for  the  1/r5  and  1/r4  asymmetries 
are  nearly  identical  until  t  ss  5.5.  In  addition,  the  decay  rate  for  the  1/r6  asymmetry  is  less  than 
that  for  the  1/r4  asymmetry  until  t  «  8.0  and  is  less  than  the  decay  rate  for  the  l/rs  asymmetry 
for  the  times  plotted.  Thus,  for  initial  conditions  tightly  confined  to  the  inner  boundary,  there  is  a 
transient  period  when  the  decay  rate  decreases  with  decreasing  disturbance  scale.  This  resembles  the 
rectilinear  problem.  The  decay  rate  eventually  reverts  back  to  increasing  with  decreasing  disturbance 
scale  though  the  cross-over  time  is  later  for  smaller-scale  disturbances. 

This  duality  may  be  qualitatively  explained  as  follows.  Initially,  most  of  the  kinetic  energy  in  the 
asymmetry  is  concentrated  near  the  inner  boundary.  The  radial  change  in  the  energy  is  so  large  that 
the  basic  state  angular  velocity  appears  constant  in  a  first  approximation.  Thus,  as  in  the  rectilinear 
problem,  decreasing  the  disturbance  scale  results  in  a  decrease  of  the  decay  rate.  However,  the  inner 
portion  of  the  asymmetry  is  symmetrized  more  rapidly  than  that  farther  out  in  the  vortex.  The 
kinetic  energy  near  the  inner  boundary  is  transferred  to  the  basic  state  much  faster  than  the  kinetic 
energy  in  the  storm’s  environment.  The  cross-over  time  occurs  when  the  inner  core  energy  decays  to 
the  same  magnitude  as  the  environmental  energy.  At  this  time,  the  basic  state  angular  velocity  no 
longer  appears  constant.  The  cross-over  time  is  later  for  smaller-scale  disturbances  because  it  takes 
longer  for  the  inner  core  kinetic  energy  to  decay  to  environmental  kinetic  energy  values. 

The  initial  vortidty  profiles  considered  thus  far  represent  forcings  that  act  nearly  continuously 
in  a  hurricane.  It  is  also  of  interest  to  investigate  transient  forcings  that  act  in  hurricanes,  such  as  an 
upper-level  short-wave  trough  or  a  tropical  easterly  wave.  As  a  proxy  for  such  forcings,  consider  a 
Gaussian  initial  vortidty  profile  of  fixed  radial  scale.  As  the  vortidty  maximum  is  moved  toward  the 
vortex  center,  the  effective  shear  and  the  decay  rate  increases  due  to  the  increase  in  the  differential 
rotation. 

The  Gaussian  initial  vortidty  profiles  considered  in  this  work  have  been  initially  upright  with 
respect  to  the  local  shear.  Consequently,  such  disturbances  have  only  had  a  decay  phase.  However, 
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migratory  short-wave  troughs,  upon  encountering  a  tropical  system,  are  generally  not  upright  with 
respect  to  the  local  shear.  As  an  example,  consider  a  trough  positioned  300  km  east  of  a  tropical 
cyclone  in  the  Northern  Hemisphere.  Let  the  trough  axis  be  oriented  30  degrees  east  of  north  such 
that  the  trough  leans  against  the  local  shear.  Recall  from  chapter  2  that  disturbances  leaning  against 
the  shear  grow  at  the  expense  of  the  basic  state.  Thus,  while  the  trough  leans  against  the  shear,  the 
trough  will  intensify  and  the  tropical  cyclone  will  weaken.  As  the  trough  grows,  it  rotates  counter 
clockwise  and  becomes  maximum  when  it  is  oriented  north-south.  Once  the  trough  is  upright  with 
respect  to  the  local  shear,  further  evolution  results  in  the  trough  weakening  and  the  tropical  cyclone 
intensifying. 

Figures  4.9, 4.10,  and  4.11  show  the  asymmetry  half-life  times  for  initial  vorticity  profiles  of  1  /r2, 
1/r3,  and  1/r8.  For  each  of  these  profiles,  the  n  =  1  asymmetry  clearly  decays  more  slowly  than  all 
higher  wavenumber  asymmetries.  The  dotted  line  in  each  plot  is  the  asymmetry  half-life  time  —  1/Scg 
for  the  given  initial  condition.  For  the  lower  wavenumbers,  —l/Seg  significantly  underestimates  the 
asymmetry  half-life  time.  For  the  1/r3  and  1/r8  profiles,  the  -1/Se«r  line  adequately  defines  the 
asymptotic  limit  of  the  half-life  time  for  the  higher  wavenumbers.  This  suggests  that  the  S„n 
definition  of  this  chapter  is  applicable  to  asymmetries  which  fall  off  in  radius  more  rapidly  than  the 
basic  state  angular  velocity.  This  conclusion  has  been  verified  for  Gaussian  initial  conditions. 

The  above  results  refute  an  assertion  made  by  CW.  Based  on  closed  form  solutions  obtained 
for  the  convection-induced  asymmetries  (cf.  CW  equations  (30)  and  (31)  ),  CW  argue  that  the 
asymmetry  damping  rate  should  be  proportional  to  the  square  of  the  wavenumber.  If  their  assertion 
was  correct,  then  as  n  — ►  oo,  the  damping  rate  should  also  become  infinite.  Thus,  the  half-life  time 
should  asymptote  to  zero.  However,  the  asymmetry  half-life  plots  show  that  for  increasing  n,  the 
time  required  for  the  energy  to  decay  to  half  its  initial  value  approaches  a  nonzero  constant  and 
does  not  decrease  as  1/n2. 

4.5  Application  to  Hurricanes 

The  definition  for  the  effective  shear  presented  in  the  previous  section  need  not  be  restricted  to 
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the  Rankine  vortex.  Here  the  effective  shear  is  evaluated  for  the  analytical  profiles  v  =  l/rQ  where 
a  =  1/3,  1/2,  and  2/3.  These  profiles  are  better  approximations  of  the  symmetric  tangential  winds 
typically  observed  in  hurricanes  than  is  the  Rankine  vortex.  The  effective  shear  is  also  calculated  'or 
the  observed  700  mb  winds  in  hurricane  Gloria  (Montgomery,  personal  communication,  1994)  which 
occurred  in  1985.  The  Gloria  data  are  nondimensionalized  using  Vm  —  57.91  ms-1  and  Rm  =  20  km 
corresponding  to  Gloria’s  maximum  tangential  wind  and  radius  of  maximum  winds,  respectively. 
The  Gloria  winds  are  given  at  10  km  intervals  and  are  considered  out  to  500  km  from  the  storm 
center.  Thus,  for  a  nondimensional  inner  boundary  set  at  a  =  1,  the  corresponding  outer  boundary 
is  at  b  =  25.  Figure  4.12  shows  the  Gloria  and  analytical  symmetric  tangential  winds. 

Since  the  effective  shear  calculation  also  requires  the  specification  of  initial  vorticity  profiles, 
three  types  of  initial  conditions  are  considered  here.  The  first  is  a  convectively  forced  asymmetry  and 
is  confined  tightly  to  the  radius  of  maximum  winds.  The  second  is  shear-induced  asymmetry  which 
is  also  confined  to  the  inner  boundary  but  extends  farther  out  into  the  vortex  than  the  convection- 
induced  asymmetry.  The  final  initial  condition  corresponds  to  a  trough  in  the  environment  of  the 
vortex  several  hundred  kilometers  from  the  vortex  center.  Mathematically,  these  initial  conditions 
are  given  by  Co  =  1/r5,  1/r2,  and  exp(-(r  -  15)2),  respectively. 

Table  4.1  shows  values  of  S*#  and  -1  /Sea  for  the  initial  voticity  profiles  and  symmetric  wind 
profiles  above.  The  values  have  been  dimensionalized  using  the  Hurricane  Gloria  data  cited  above. 
The  asymmetry  half-life  times  given  by  -1/Seff  are  underestimates  of  the  decay  rate.  Examination 
of  figures  4.9,  4.10,  and  4.11  indicate  the  n  =  1  asymmetry  decays  approximately  seven  times  slo  3r 
than  the  limiting  time-scale.  Assuming  that  the  n  =  1  decay  rate  approximates  the  decay  rate  for  an 
atmospheric  disturbance,  the  asymmetries  forced  in  the  inner  region  of  a  hurricane  give  up  half  their 
initial  energy  in  35-70  minutes  for  the  analytical  symmetric  winds.  The  table  suggests  that  similar 
asymmetries  in  Hurricane  Gloria  gave  up  half  their  initial  energy  in  40-60  minutes.  The  Gaussian 
initial  vorticity  here  represents  a  trough  300  km  from  the  storm  center.  This  feature  decays  to  half 
its  initial  energy  in  17-35  hours  for  the  analytical  symmetric  winds.  The  table  suggests  that  similar 
asymmetries  in  Hurricane  Gloria  gave  up  half  their  initial  energy  in  14  hours. 
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4.6  Summary 


Wave  disturbances  on  a  symmetric  vortex  have  been  shown  to  favor  the  lower  wavenumber 
asymmetries.  Disturbances  evolve  in  a  similar  fashion  to  their  rectilinear  counterparts  for  changes 
in  the  wavenumber  and  in  time.  However,  the  disturbance  evolution  for  changing  radial  structure  is 
complicated  by  the  change  in  the  basic  state  angular  velocity  with  radius.  In  addition  to  the  radial 
scale,  the  radial  location  of  the  initial  asymmetry  must  be  known. 

For  asymmetries  forced  in  the  inner  core  region  of  the  vortex,  decreasing  the  radial  scale  of  the 
asymmetry  increases  the  effective  shear  and  the  decay  rate  when  the  environmental  portion  of  the 
kinetic  energy  is  comparable  to  the  inner  core  portion.  When  the  asymmetry  is  tightly  confined  to 
the  inner  boundary,  the  kinetic  energy  near  the  inner  core  dominates.  The  effective  shear  appears 
constant  and  further  narrowing  of  the  radial  scale  results  in  a  decreasing  decay  rate.  This  transient 
behavior  lasts  until  the  kinetic  energy  near  the  inner  core  is  reduced  to  levels  comparable  to  the 
environmental  kinetic  energy.  For  asymmetries  forced  in  the  environment  of  the  vortex,  increasing 
the  radial  scale  of  the  asymmetry  or  moving  a  fixed  scale  asymmetry  towards  the  inner  boundary 
results  in  an  increase  of  the  effective  shear  and  the  decay  rate. 

Since  the  shear  varies  with  radius,  an  integrated  measure  of  the  shear  is  needed  to  estimate  the 
global  evolution  of  asymmetries.  The  shear  definition  presented  here  considers  both  the  basic  state 
angular  velocity  and  the  initial  vorticity  profile.  Noting  that  the  reciprocal  of  the  rectilinear  shear 
defined  the  limiting  time  for  a  perturbation  to  decay  to  half  its  initial  energy,  the  reciprocal  of  the 
effective  shear  was  expected  to  define  the  limiting  asymmetry  half-life  time.  This  time  scale  signifi¬ 
cantly  underestimated  the  decay  rate  for  the  lower  azimuthal  wavenumbers  for  all  initial  conditions 
considered.  However,  the  decay  time  scale  -1  /Sea  reasonably  represented  the  asymptotic  limit  of 
the  asymmetry  half-life  times  for  the  higher  wavenumbers  as  long  as  the  initial  condition  fell  off  in 
radius  more  rapidly  than  the  basic  state  angular  velocity. 
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Table  4-1:  Effective  shear  values  and  limiting  asymmetry  half-life  times  for  various  symmetric  wind 
profiles.  Symmetric  wind  data  from  the  700  mb  level  of  Hurricane  Gloria  are  also  considered.  All 
values  have  been  dimensionalized  using  Gloria’s  maximum  azimuthal  wind  (57.91  ms-1)  and  the 
radius  of  maximum  winds  (20  km). 

Co 


1/r5 

1/r2 

g(r— 15)a/10 

l/rX/3 

-2.81  x  10~3  s"1 

-1.77  x  10~3  s"1 

— 1.09  x  10“4  s- 

5.94  min 

9.44  min 

2.56  hr 

1/r1/2 

-3.06  x  10~3  s"1 

-1.86  x  10~3  s"1 

-7.83  x  10-®  s— ' 

5.45  min 

8.96  min 

3.55  hr 

1/r2/3 

-3.29  x  10~3  s"1 

-1.94  x  10~3  8_1 

-5.59  x  10~5  s- 

5.07  min 

8.58  min 

4.97  hr 

Gloria 

-2.89  x  10-3  s"1 

-1.81  x  10~3  s"1 

-1.29  x  10~4  s-! 

5.76  min 

9.23  min 

2.15  hr 
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Normalized  Kinetic  Energy 
Bounded  Vortex:  a=l,  b=10 


TIME  (HRS) 


Fig.  4.2:  Figure  4.(a)  from  the  Carr  and  Williams  (1989)  paper.  It  is  a  plot  of  the  normalized 
kinetic  energy  as  a  function  of  time  for  the  four  initial  conditions  developed  in  their  paper.  The 
curves  are  the  energy  decay  for:  convection-induced  (n  =  2,  dot;  n  =  1  dash),  motion-induced 
(chaindot),  and  shear-induced  (chaindash)  asymmetries. 


Time  (HRS) 


Fig.  4.3:  Normalized  kinetic  energy  as  a  function  of  time  for  the  four  initial  conditions  developed 
by  Carr  and  Williams  (1989).  The  curves  are  the  energy  decay  for:  convection-induced  (n  =  2,  dot; 
n  =  1,  dash),  motion-induced  (dot  dash),  and  shear-induced  (dot  dot  dot  dash)  asymmetries. 


Perturbation  Streamfunction 
Bounded  Vortex:  a=l,  b=10 
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Fig.  4.4  :  Perturbation  streamfunction  field  for  Co  =  1/r3  for  the  bounded  Rankine  vortex.  The 
columns  show  contours  at  times  t  =  0.0,  t  =  3.6,  and  t  =  7.2,  respectively.  Panels  (a)-(c)  show 
contours  for  n  =  1  where  the  <  ontour  interval  is  1.10  x  10-2.  Panels  (d)-(f)  show  contours  for  n  =  2 
where  the  contour  interval  is  5.67  x  10-3.  Panels  (g)-(i)  show  contours  for  n  =  3  where  the  contour 
interval  is  3.31  x  10“3. 
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Fig.  4.5:  Perturbation  vortidty  field  corresponding  to  figure  4.4.  Note  the  different  axis 
The  contour  interval  is  6.90  x  10-2. 
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Normalized  Kinetic  Energy 
Bounded  Rectilinear  Channel:  a=l,  b=10 


Fig.  4.6:  Normalized  kinetic  energy  for  the  rectilinear  problem  with  =  — 1/(3/  -  ll)3.  The  curves 
are  the  energy  decay  for:  k  =  1  (dot),  k  =  2  (dash),  k  =  3  (dot  dash).  The  solid  curve  is  the  energy 
decay  for  the  upright  plane  wave  given  by  1/(1  + 12). 
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Normalized  Kinetic  Energy 
Bounded  Vortex:  a=l,  b=10 


Fig.  4.7:  Normalized  kinetic  energy  for  the  bounded  Rankine  vortex  with  Co  =  1/r3.  The  curves 
are  the  energy  decay  for:  n  =  1  (dot),  n  =  2  (dash),  n  =  3  (dot  dash).  The  solid  curve  is  the  limiting 
energy  decay  given  by  1/[1  +  (Segt)3]  where  Stg  =  -1.01. 


67 


Fig.  4.8:  Normalized  kinetic  energy  for  the  bounded  R&nkine  vortex  for  n  =  1.  The  curves  are  the 
energy  decay  for:  Co  =  1  /r3  (dot),  Co  =  1/r4  (dash),  Co  =  1  /r6  (dot  dash),  and  Co  =  1/r6  (dot  dot 
dot  dash). 
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Asymmetry  Half-Life 
Bounded  Vortex:  a=l,  b=10 


I 

j 

i 

I 


Fig.  4.9:  Time  required  for  an  asymmetry  to  decay  to  half  its  initial  energy  as  a  function  of 
wavenumber  for  Co  =  1/r2.  The  dotted  line  is  the  limiting  case  given  by  -1/S«ff  where  Seff  =  -0.74. 
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Asymmetry  Half-Life  Time 


Asymmetry  Half-Life 
Bounded  Vortex:  a=l,  b=10 


N 


Fig.  4.10:  Same  as  figure  4.9  except  Co  =  1  /r3  and  Seg  =  -1.01. 


Asymmetry  Half-Life 
Bounded  Vortex:  a=l,  b=10 
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Fig.  4.11:  Same  as  figure  4.9  except  Co  =  1/r5  and  Seff  =  -1.34. 
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CHAPTER  5 


CONCLUSION 


5.1  Summary 

The  inviscid  nature  of  disturbance  evolution  in  sheared  flows  has  been  investigated.  The  problem 
has  been  studied  within  the  framework  of  nondivergent  vorticity  dynamics.  Particular  emphasis  has 
been  placed  on  understanding  the  wavenumber  and  spatial  structure  dependencies  in  perturbation 
evolution. 

A  truly  inviscid  mechanism  which  favors  the  decay  of  high  wavenumber  perturbations  over 
low  wavenumber  perturbations  was  identified.  Further,  the  wavenumber  one  perturbation  decayed 
significantly  slower  than  all  other  perturbations.  In  the  development  of  AB  theory,  SM  appealed  to 
both  observational  and  computational  evidence  of  a  wavenumber  selection  mechanism  as  the  basis 
for  the  theory.  The  results  presented  here  further  elucidate  the  wavenumber  selection  mechanism 
and  provide  further  theoretical  support  for  the  validity  of  AB  theory. 

The  spatial  structure  of  the  initial  conditions  was  shown  to  play  a  significant  role  in  disturbance 
evolution.  The  spatial  scale  controlled  how  rapidly  a  disturbance  sheared  and  subsequently  decayed. 
In  the  swirling  problem,  the  radial  location  of  the  initial  condition  was  an  important  consideration 
as  the  shear  decreased  from  the  radius  of  maximum  winds.  These  factors  led  to  the  definition  of 
an  effective  shear  which  accounted  for  the  interaction  of  the  symmetric  basic  state  and  the  initial 
vorticity  profile.  As  defined  here,  the  effective  shear  was  shown  to  be  an  adequate  estimate  for  the 
higher  wavenumbers  as  long  as  the  radial  variation  of  the  initial  condtion  was  greater  than  that  of 
the  symmetric  basic  state.  However,  the  effective  shear  was  a  significant  underestimate  for  the  lower 
wavenumbers. 

Using  the  effective  shear,  decay  half-life  times  were  estimated  on  a  wavenumber  by  wavenumber 
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basis.  The  results  clearly  demonstrated  that  the  decay  half-life  times  do  not  decrease  as  1  /n2  which 
refuted  the  CW  assertion  that  the  damping  rate  is  proportional  to  the  square  of  the  azimuthal 
wavenumber.  The  effective  half-life  times  for  wavenumber  one  were  found  to  be  approximately 
seven  times  slower  than  those  forecast  by  the  effective  shear.  This  was  used  to  estimate  half-life 
times  for  tangential  wind  profiles  that  are  more  representative  of  hurricanes  outside  their  radius 
of  maximum  winds.  These  half-life  times  seemed  to  be  reasonable  estimates  for  the  asymmetric 
forcings  simulated. 

5.2  Suggested  Further  Work 

The  wind  profiles  considered  did  not  have  a  basic  state  vorticity  gradient  and,  thus,  excluded 
discrete  normal  modes  from  the  solution  set.  While  the  interaction  of  the  continuous  spectrum 
and  the  discrete  normal  modes  has  been  considered  for  the  Eady  model  (Pedlosky,  1964;  Farrell, 
1984),  this  topic  requires  further  investigation  in  rapidly  rotating  vortices.  An  analytical  method 
that  extends  the  Rankine  vortex  into  a  uniformly  rotating  core  is  briefly  considered  in  Appendix  A. 
In  addition,  a  basic  state  vorticity  gradient  may  be  viewed  as  an  effective  vortex  0.  CW  suggested 
that  this  effective  0  would  introduce  a  retrogression  that  could  counteract  the  rate  of  perturbation 
tilting  induced  by  the  differential  rotation.  Further  work  is  needed  to  examine  this  hypothesis. 

Sutyrin  characterized  the  transfer  of  energy  from  the  azimuthal  perturbations  to  the  circularly 
symmetric  basic  state  as  being  analogous  to  the  energy  cascade  to  larger  scales  in  two-dimensional 
turbulence.  Moreover,  he  associated  the  growth  of  the  basic  state  potential  vorticity  gradient  with 
the  transfer  of  enstrophy  to  smaller  scales.  These  ideas  merit  further  investigation  and  quantification. 
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APPENDIX  A 


AXISYMMETRIZATION  IN  UNBOUNDED  RANKINE  FLOW 


In  this  appendix,  the  asymmetric  dynamics  of  an  unbounded  Rankine  vortex  are  briefly  consid¬ 
ered.  The  nondimensional  linearized  vorticity  equation  from  chapter  4  is  again  used  to  represent  the 
system.  The  inner  boundary  is  now  at  the  origin  while  the  outer  boundary  is  removed  to  infinity. 
It  is  assumed  that  all  quantities  are  bounded  at  the  origin  and  vanish  at  infinity. 

When  including  the  region  inside  the  radius  of  maximum  winds,  the  Rankine  profile  is  given  by 


{r/o,  r  <  n; 
a/r,  r  >  o; 


(A.1) 


where  a  is  the  nondimensional  radius  of  maximum  winds  and  r  is  the  nondimensional  radius.  The 
corresponding  basic  state  vorticity  profile  is 


C  =  [  2/°’  r  <  °; 
s  \  0,  r  >  o; 

while  the  absolute  vorticity  profile  is  given  by 


(A.  2) 


The  quantity  R~l  is  the  inverse  Rossby  number  defined  by  fZ"1  =  fRm/Vm  where  /  is  the  Coriolis 
parameter  and  Vm  is  the  symmetric  wind  at  the  radius  of  maximum  winds  (fin).  Figure  A.l  shows 
S  and  i}  for  the  full  Rankine  vortex. 

The  discontinuity  in  the  mean  state  vorticity  at  r  =  a  effectively  introduces  another  boundary 
to  the  system.  Since  this  boundary  is  interior  to  the  fluid,  the  kinematic  and  dynamic  boundary 
conditions  must  be  satisfied  at  the  disturbed  interface  r  =  a+e,  where  e  is  the  interface  displacement. 
The  kinematic  boundary  condition  requires  that  the  normal  velocity  be  continuous  while  the  dynamic 
boundary  condition  requires  that  the  pressure  be  continuous  at  r  =  a  +  e,  respectively.  Consistent 
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with  the  linearization  employed  below  the  matching  conditions  will  be  evaluated  at  r  =  a.  The 
evolution  of  the  disturbed  interface  may  still  be  determined  by  integrating 


(d_  C6_\ 

\dt  +  rd\) 


(  —  u 


M-4) 


at  r  =  a  once  tt'  is  found. 

For  the  full  Rankine  profile  (A.2)  the  linearized  vorticity  equation  of  chapter  4  is  modified  to 


(d_  £^_\ 

+  r  d\) 


f 


=  0, 


r  a. 


(A.  5) 


To  solve  (.4.5),  the  discontinuity  in  the  basic  state  vorticity  at  r  =  a  must  be  accounted  for.  Since 
the  problem  is  linear,  the  superposition  principle  may  be  used  to  decompose  the  solution  into  two 
parts  by  letting  C  —  C*  +  Ci  where  £,  is  defined  to  be  smooth  for  all  r  and  £i  accounts  for  the 
discontinuity  in  the  basic  state  vorticity  at  r  =  a.  The  vorticity  equation  (.4.5)  is  then  split  into  two 
parts 


<A6i> 


Equation  (.4.6a)  is  formally  identical  to  the  system  solved  in  chapter  4,  but  with  the  boundary 
conditions  cited  above.  The  corresponding  solution  in  Fourier  space  is 

OO 

Mr,  t)  =  j  G(r,  pX.o(p)e-'nBt/p  pdp  (41.7) 

o 

where  the  appropriate  Green’s  function  is 


G(rip)  =  _J_/p~n’-">  0<r<p 
V  2n  \  pnr~n,  p<r<  oo, 


and  C.o(p)  is  the  smooth  component  of  C  at  time  t  =  0. 
The  Fourier  space  equivalent  to  (AM)  is 


f  d  tnf>\  ; 

U+-)<‘=0' 


r  £  a. 


(A.&) 


(A.9) 
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Anticipating  that  the  solution  to  (A. 9)  will  yield  the  discrete  normal  modes  which  are  irrotational 
on  both  sides  of  the  discontinuity,  (i  is  assumed  separable  of  the  form 

Ci  =  l(t)6{r  -  a).  (A. 10) 

Here,  7  is  an  undetermined  temporal  multiplier  for  Ci  and  6(r  -  a)  is  the  Dirac  delta  function.  In 
terms  of  the  streamfunction,  equation  (A.10)  is  written 


VJV»i  =  7 (t)6(r  -  a). 


(A.ll) 


The  streamfunction  is  also  assumed  separable  of  the  form  rpi  =  7(t)4'i(r).  Thus,  equation  (A.ll) 
becomes 

For  r  ^  o,  equation  (A.12)  is  Euler’s  equation.  Two  conditions  are  needed  to  match  the  solutions 
in  each  region  across  r  =  a.  The  first  is  the  kinematic  boundary  condition  requiring  that  the  radial 
velocity,  u,  be  continuous  at  r  —  a.  Thus,  the  Fourier  streamfunction  amplitude  must  be  continuous 
across  r  =  a.  The  second  condition  results  from  integrating  (A.12)  over  a  small  interval  that  includes 
r  =  a.  This  yields  the  following  jump  condition  for  $1 


3$i  + 

■5T(*  >  -  ~5T(o  >  =  L 


(A-13) 


Applying  the  boundary  conditions,  and  the  continuity  and  jump  conditions  at  r  =  a,  yields 

A.  - _ 5L  /  °"nr">  0  <  r  <  o  ,A  14v 

Wl  2n  \  onr_n,  a  <  r  <  00.  ' 

To  complete  the  derivation,  7  must  be  determined.  The  remaining  constraint  is  the  dynamic 
boundary  condition  which  requires  that  the  pressure  be  continuous  at  r  =  a.  In  Fourier  space,  the 
azimuthal  momentum  equation  is  given  by 

«ii  <n 

(A-15) 


3ti  in€  „  in. 

■57  +  —  v  +  fpi  =  p. 
eft  r  r 


Evaluating  (A.  15)  on  each  side  of  r  =  a  and  subtracting  gives 


j|[t>(a+)  -  v(o")j  +  ^(6(o+)  -  t)(a  )]  -  ^u(o)  =  0. 
eft  a  a 


(A.16) 
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In  terms  of  the  total  streamfunction,  equation  (A.  16)  becomes 


d_ 

dt 


in 

+  — 
a 


<%>,  +x  dip  _ 
* <a  > -  *  <a  > 


+  ^(a)  =  0.  (All) 


Now,  from  the  superposition  principle,  t£>  =  ip,  +  7^1 ,  where  ip,  and  its  derivatives  are  everywhere 
smooth  by  construction.  Since  is  continuous  but  has  a  unit  jump  in  its  derivative  across  r  =  o, 
equation  (A.17)  simplifies  to 

§  +  ^(n-l)7  =  -~ip,(a,t),  (A.18) 

a  first-order  linear  differential  equation  for  7.  Upon  multiplying  through  by  the  integrating  factor 
exp[t(n  -  l)t/o]  and  substituting  for  ip,(a,t),  equation  (A.18)  becomes 

OO 

|  ^  J  G(a ,  pJC.oeK"-1  >/-**/*  ^  (A.  19) 

0 

Integrating  in  time  and  then  multiplying  through  by  exp[— t(n  —  l)t/a]  gives 

0 

where  C\  is  the  constant  of  integration  at  t  =  0.  Equation  (A.20)  then  yields 

’  Pdp 


=  f  G(tt»P)C«o  -inW/t 

7W  a  J  (n  —  1  —  ana/p) 


2n  7  G(a,pK,0 

0  a  J  (n  —  1  -  ana  jp) 
0 


c-*(n-l)t/a 


(A.21) 


where  70  is  the  initial  amplitude  of  the  normal-mode  (Rossby)  wave.  The  Fourier  streamfunction 
amplitude,  ip,  is  thus 


OO 

4  =  J  G(r,p)Ue-in9t/p  pdp 

0 

+  *■ 


-*1 


■  ana /p) 
G(a,p)  C» 


«%  ,  f  G(a,  p)C»o  . 
2n  *  J  (n  -  1  -  ana  Ip) 


e-i(n~l)t/a 


(A.  22) 


where  G  is  given  by  (A.8)  and  4ri  is  given  by  (A. 14).  To  obtain  the  physical  space  streamfunction, 
the  inverse  Fourier  transform  must  be  applied  to  (A.22). 
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Equation  ( A.22 )  is  interpreted  as  follows.  The  first  line  may  be  tentatively  identified  with  the 
continuous  spectrum  solution  and  is  the  analogue  of  the  solution  presented  in  chapter  4.  The  second 
line  is  a  conversion  term  that  transfers  some  of  the  energy  in  the  continuous  spectrum  solution  into 
the  discrete  mode,  which  is  the  third  line  of  the  equation.  Interestingly,  even  with  no  normal  mode 
component  initially  (7<>  =  0),  the  continuous  spectrum  solution  will  always  project  onto  the  normal 
mode  at  later  times  (cf.  Farrell,  1982).  Note  that  for  n  =  1,  the  discrete  mode  does  not  rotate  on 
the  vortex,  but  rather  represents  the  translation  of  the  basic  state  vortex.  For  n  ^  1,  the  discrete 
modes  rotate  slower  than  the  vortex  and  represent  retrogressing  Rossby  waves  at  r  =  a  +  e. 

Edwards  (1994)  demonstrated  that  this  model  can  be  extended  to  an  arbitrary  number  of 
regions  of  constant  basic  state  vortidty.  Sutyrin  (1989)  developed  the  quasi-geostrophic  shallow 
water  equivalent  of  the  multi-region  model.  However,  his  formulation  did  not  explicitly  describe  the 
interaction  between  the  continuous  and  discrete  spectrum  solutions.  Extension  of  the  above  model 
to  a  quasi-geostrophic  shallow  water  system  is  of  interest  as  it  may  prove  useful  in  idealized  studies 
of  hurricane  track. 
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APPENDIX  B 


MODEL  VERIFICATION 


In  this  appendix,  the  verification  of  the  numerical  model  is  briefly  considered.  Since  the  merid¬ 
ional  or  radial  disturbance  scale  generally  changes  with  time,  an  estimate  of  when  the  disturbance 
scale  becomes  smaller  than  the  model  resolution  is  desired.  In  the  rectilinear  problem  of  chapter  2, 
the  meridional  wavenumber  of  an  upright  plane  wave  perturbation  was  found  to  be  l  —  —kSt.  Thus, 
a  meridional  wavelength  was  given  by  2x/kSt.  Assuming  the  model  resolves  perturbations  that  have 
scales  greater  than  or  equal  to  twice  the  grid  spacing,  the  model  no  longer  resolves  the  perturbations 
when 

(PA) 


2x 

kS,*2A'1 


where  Ay  is  the  grid  spacing,  k  is  the  zonal  wavenumber,  and  5  is  the  shear.  Thus, 

ic 


St> 


kAy 


(B- 2) 


gives  an  estimate  of  when  the  perturbations  are  no  longer  resolvable.  In  chapter  2  and  3,  the  worst 
model  resolution  was  Ay  =  0.056,  corresponding  to  361  grid  points  and  integration  limits  at  ±10. 
For  k  —  3,  the  model  resolution  starts  becoming  inadequate  for  St «  20. 

For  the  swirling  problem,  SM  have  shown  that  the  change  in  the  radial  wavenumber  with  time 
is 

dk  dft  ,DO, 

*  “  -*•  (B3) 

Here,  k  is  the  local  radial  wavenumber,  n  is  the  azimuthal  wavenumber,  and  ft  is  the  basic  state 
angular  velocity.  For  the  Ranldne  profile  (A.  1),  equation  (B. 3)  is 

§—£•  (*-4> 

Thus,  for  a  locally  upright  asymmetry 

(B.  5) 


,  2  nt 

*=—■ 
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The  model  resolution,  then,  starts  becoming  inadequate  when 


27T 

—  <  2Ar 
k 


(B.  6) 


or,  when 


t  > 


*r3 


(B.7) 


~  2nA r 

In  chapter  4,  the  worst  model  resolution  was  Ar  =  0.05,  corresponding  to  boundaries  a  —  1  and 
6  =  10  with  181  grid  points.  Thus,  at  the  radius  of  maximum  winds  for  n  —  3,  the  asymmetries 
start  becoming  smaller  than  the  model  resolution  when  t  as  10.  However,  at  r  =  5  and  n  =  3,  the 
asymmetries  start  becoming  unresolvable  when  t  as  1300. 

To  be  consistent  with  the  analytical  system,  the  model  results  must  satisfy  the  momentum 
equations  to  within  discretization  error.  Since  the  trapezoidal  rule  was  used  in  this  work,  the 
discretization  error  should  be  proportional  to  the  model  resolution  squared,  i.e.,  model  output  must 
be  at  least  quadratically  convergent.  The  procedure  used  to  evaluate  the  model  output  for  the 
swirling  problem  is  briefly  presented  below. 

For  the  bounded  Rankine  vortex  (4.18),  model  verification  was  performed  with  the  Fourier  space 
representation  of  the  momentum  equations  (4.8a)  and  (4.86) 


du  in . 

a.*?-*’- 

dv  in. 

-  +  -v  +  1u  = 


dp 
dr' 
in . 

-Tp. 


(B.&a) 
{B.  86) 


From  (£.86),  p  was  found  for  all  radial  values  at  a  fixed  time.  Using  centered  difference  approxima¬ 
tions  for  all  differentiated  quantities,  the  radial  gradient  of  p  was  calculated  and  compared  to  the 
left  side  of  (B.&a).  For  the  shear-induced  asymmetry  from  CW  (£„  =  — i/r2,  n  =  2)  and  model 
resolutions  of  Ar  =  At  =  5  x  10~2  at  t  =  9,  the  difference  between  the  left  and  right  sides  of  (B.&a) 
at  r  =  2.3  was  -2.54  x  10~s.  For  double  and  quadruple  the  above  model  resolution,  the  difference 
was  —6.34  x  10~6  and  -1.59  x  10~8,  respectively.  In  each  case,  as  the  model  resolution  doubled, 
the  difference  decreased  by  a  factor  of  4.  For  the  same  model  resolutions  at  r  =  1.7  and  t  =  5.9, 
the  differences  were  9.02  x  10~®,  2.01  x  10-9,  and  4.87  x  10~7,  respectively.  The  model  output  was 
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checked  for  a  variety  of  times,  radial  locations,  and  initial  vorticity  profiles.  In  all  cases,  the  model 
output  was  found  to  be  quadratically  convergent. 
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